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This  summary  report  is  a  collection  of  four  separate 
progress  reports  prepared  under  three  contracts,  which  are 
all  sponsored  by  the  Office  of  Naval  Research  in  Arlington, 
Virginia.  This  report  contains  the  results  of  investiga¬ 
tions  into  the  app.  ation  of  the  maximum  entropy  method 
(MEM),  a  high  resolution,  frequency  and  wavenumber  estima¬ 
tion  technique.  The  report  also  contains  a  description  of 
two,  new,  stable  ,  high  resolution  spectral  estimation  tech¬ 
niques  that  is  provided  in  the  final  report  section.  Many 
examples  of  wavenumber  spectral  patterns  for  all  investiga¬ 
ted  techniques  are  included  throughout  the  report. 

The  maximum  entropy  method  is  also  known  as  the  maxi¬ 
mum  entropy  spectral  analysis  (MESA)  technique,  and  both 
names  are  used  in  the  report.  Many  MEM  wavenumber  spectral 
patterns  are  demonstrated  using  both  simulated  and  measured 
radar  signal  and  noise  data.  Methods  for  obtaining  stable 
MEM  wavenumber  spectra  are  discussed,  broadband  signal  de¬ 
tection  using  the  MEM  prediction  error  transform  (PET)  is 
discussed,  and  doppler  radar  narrowband  signal  detection  is 
demonstrated  using  the  MEM  technique.  It  is  also  shown 
that  MEM  cannot  be  applied  to  randomly  sampled  data. 

The  two  new,  stable,  high  resolution,  spectral  esti¬ 
mation  techniques  discussed  in  the  final  report  section, 
are  name'  the  Wiener -King  and  the  Fourier  spectral  estimation 
techniques.  The  two  new  techniques  have  a  similar  deriva¬ 
tion  based  upon  the  Wiener  prediction  filter,  but  the  two 
techniques  are  otherwise  quite  different.  Further  develop¬ 
ment  of  the  techniques  and  measurement  of  the  technique  spec¬ 
tral  characteristics  is  recommended  for  subsequent  investiga¬ 
tion. 
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STABLE  MESA  ANTENNA  PATTERNS 


I.  INTRODUCTION 

A.  Wavenumber  Power  Spectra 

Power  spectral  analysis  techniques  are  applicable  to  the 
processing  of  spatial  multi-channel  antenna  data,  since  computed 
antenna  patterns  are  actually  wavenumber  power  spectra,  P  (k )  .  The 
wavenumber  Ms  a  function  of  9,  the  signal  angle  of  incidence  to 
the  antenna  as  follows: 

k  =  (2tt/X)  SIN  (9) 

Because  an  antenna  array  is  a  collector  of  spatially  sampled  data  , 
any  power  spectral  technique  which  is  designed  for  such  discrete 
data  sets  is  applicable  for  computing  antenna  patterns. 

B.  High  Resolution  Power  Spectra  1  Techniques 

During  recent  years  several  high  resolution  power  spectral 
techniques  have  been  developed  (or  rediscovered)  for  use  with  dis¬ 
crete  data  sets.  Some  of  the  techniques  are:  the  maximum  entropy 
method  (1),  the  autoregressive  model  (2),  the  nv  ving  average  model  (3)  , 
the  Yule-Walker  technique  (4),  and  the  maximum  riklihood  method  (5). 

Of  these  techniques  the  maximum  entropy,  autoregressive,  and  Yule-Walker 
techniques  are  all  pole  models,  the  moving  average  technique  is  an  all 
zero  model,  and  maxiumum  liklihood  is  only  a  criteria  function  applicable 
to  any  model.  Most  of  these  methods  are  described  in  a  tutoral  review 
article  (6).  These  particular  methods  have  also  been  investigated 
and  compared  in  two  reports  (7) ,  (8)  in  which  the  best  results  were 

achieved  with  the  maximum  entropy  method. 

More  conventional  high  resolution  Fourier  methods  have 
also  been  recently  developed  (9) ,  (10) ,  (11) ,  but  have  not  been  so 

thoroughly  investigated.  As  a  consequence ,  in  this  paper  several 
methods  are  investigated  for  applying  the  maximum  entropy  spectral 
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analysis  technique  (MESA)  to  the  processing  of  spatial,  uniformly 
sampled  data. 

C .  Application  of  MESA 

Several  methods  of  applying  MESA  to  spatial  data  are  inves¬ 
tigated.  In  particular  Spatial  daca,  which  is  simulated  for  an  8 
element  linear  antenna  array,  is  processed  with  MESA  and  the  Burg 
technique  (1).  The  Burg  technique  j.s  a  recursive  method  for  eval¬ 
uating  the  MESA  filter  weights,  which  substantially  reduces  the 
number  of  calculations  required  of  the  more  conventional  inverse 
matrix  evaluation  method. 

MESA  antenna  patterns  (wavenumber  spectra)  may  be  computed 
upon  the  collection  of  the  set  of  8  spatial  data  samples  at  any 
instant  of  time.  Such  "snapshot"  patterns  are  inherently  inconsistent 
and  unstable.  However  it  is  possible  to  compute  stable  MESA  antenna 
patterns  using  one  of  several  stablization  techniques.  It  remains 
only  to  determine  which  technique  provides  sufficient  stability  for 
an  acceptable  averaging  period  without  destroying  the  desired  high 
resolution  property  which  is  characteristic  of  MESA.  Averaging 
techniques  which  are  investigated  employ  a  time  average  of  one  of 
the  follwoing  sets  of  variables: 

a. )  filter  weights 

b. )  prediction  errors 

c. )  covariance  matrix 

d. )  "snapshot"  patterns 

And  as  an  alternative  to  averaging,  stablizatio  may  also  be  achieved 
with  use  of  time  adaptive  filter  weights.  In  particular  a  set  of 
adaptive  filter  weights,  which  are  defined  as  proportional  to  the 
prediction  error  (12),  are  utilized  in  conjunction  with  a  proportion¬ 
ality  constant  (convergence  parameter)  tc  comprise  a  stable,  adaptive 
MESA  processing  technique. 

II.  THE  MAXIMUM  ENTFOPY  METHOD 

"'he  MESA  technique,  as  the  name  implies ,  originated  (1)  by 


maximizing  the  entropy  of  a  signal  mixed  with  noise.  However,  the 
same  filter  weights  may  also  be  derived  (13)  by  whitening  the  Weiner 
prediction  error  filter  as  specified  for  discreet  data  samples  (14). 
The  resulting  maximum  entropy  wavenumber  spectra  P(k)  is  given  as 
follows : 


P(k) 
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N 


N 

1  +  L  y 
n~l 


N 
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exp (iknAx) 
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(1) 


where 

N  =  number  of  filter  weights  (1-N-M) 

M  =  number  of  data  samples 
PN  =  total  noise  power 

AX  =  antenna  element  spacing 
YN 

n  =  nth  prediction  error  filter  weight  of  a  set  of  N  weights 

The  variables  of  eqn.  (1),  which  are  computed  using  a 
set  of  equations  known  as  the  "Burg  technique",  are  listed  as  follows: 


a . 


Total  Noise  Power  Pv 


Pi  =  r-  (r_  is  the  data  set  autocorrelation  function) 
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where  the  Kth  data  sample  =  X 


for  ( 1-K-M ) 


III.  MESA  SNAPSHOTS 


It  .is  possible  to  compute  an  antenna  pattern  with  MESA 
using  only  one  set  of  M  data  samples  all  recorded  at  the  same  instant 
of  time.  For  example,  consider  one  set  of  8  data  samples  collected 
with  8  unformly  spaced  antennas.  The  Burg  technique  equations,  eqns. 
(1-6)  are  initially  evaluat jd  for  N=1  and  M=8 ,  and  then  evaluated 
repeatedly  for  increasing  unit  incremental  values  of  N  up  to  the 
desired  value  of  M  provided  that  (lfN£7), 
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However  the  final  value  of  N  must  be  such  that  (NSiM)  where  NS  is 
the  number  of  signals  present  in  the  given  data  set  X„. 

Consider  a  MESA  snapshot  pattern  evaluated  for  N=4,M=8 
where  the  Burg  technique  equations  are  evaluated  repeatedly  for 
(1SNS4) .  A  MESA  snapshot  pattern  of  one  signal  incident  at  +10 
degrees  'o  degrees  is  broadside  to  the  antenna)  and  a  signal-to- 
noise  ratio  of  15  dB  is  shown  in  Fig.  la.  The  8  data  points  contain 
Guassian,  white  noise,  simulated  using  a  set  of  8  random  numbers 
computed  for  a  generator  "seed”  value  of  1  (IR=1).  Another  MESA 
snapshot  shown  in  Fig.  lb,  is  computed  using  a  different  set  of  8 
random  numbers  for  which  IR=2 . 

The  single  signal  is  located  accurately  at  +10  degrees 
(within  -0.5  degrees)  in  both  MESA  snapshots  of  Fig.  1.  The  side 
peaks,  which  are  randomly  located,  occur  at  different  positions 
in  the  two  snapshots.  The  total  number  of  peaks,  which  represent 
the  poles  of  eqn.  (1) ,  is  always  less  or  equal  to  the  value  of  N, 
the  number  of  filter  weights.  Since  the  two  independent  data  sets 
usedrin  the  computed  antenna  patterns  of  Fig.  1  are  considered  to  be 
recorded  at  two  different  instances  of  time,  MESA  snapshots  are 
clearly  time  variant  when  computed  with  short  (M=8)  data  sets.  It 
is  evident  that  some  stablizing  technique  is  needed  in  the  application 
of  MESA  to  nhOit  data  sets,  so  that  computed  .MESA  antenna  patterns 
are  invariant  and  repeatable  in  time  for  stationary  data. 

IV.  MESA  INSTABILITIES 

Besides  the  side  peak  location  instability  depicted  by 
Fig.  1,  another  instability  associated  with  MESA  is  the  inaccurate 
representation  of  signal  peaks.  Signal  peaks  may  not  be  accurately 
locate d  at  very  Low  signal- to-noise  ratios  or  when  other  signals 
are  present  at  adjacent  angles.  Nearby  signals  cause  distortion 
in  both  signal  location  and  in  relative  signal  peak  height.  It 
has  oeen  noted  (+5)  that  isolated  MESA  signal  peaks  are  not  linearly 
related  to  the  signal-tc-noise  ratios,  but  the  actual  relationship 
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Fig.  lb  -  MESA  snapshot,  one  signal,  IR  2 


has  not  been  demonstrated.  Both  the  power  lever  and  relative  phase 
of  adjacent  signals  affect  the  accuracy  of  the  MESA  power  spectra. 
Such  signal  instabilities  have  been  noted  previously  (14). 

The  noise  field  associated  with  a  multi-channel  antenna 
array  may  also  be  a  source  of  signal  distortion.  Split  signal  peaks 
are  a  common  problem  with  MESA  snapshots  (16).  An  example  of  how 
such  problems  arise  is  depicted  in  Figs.  2  and  3.  The  MESA  snapshot 
of  Fig.  2  is  computed  using  only  3  complex  noise  data  points  (no 
signal  present)  and  a  5  point  filter  (N=5).  The  MESA  snapshot  of 
Fig.  2  has  two  large  noise  peaks  near  0  and  10  degrees.  When  a 
signal  having  a  10  dB  signal-to-noise  ratio  is  introduced  into  the 
same  noise  field  (IR=5),  at  an  angle  of  +5  degrees,  the  signal  is 
ambiguously  represented  in  the  MESA  snapshot  of  Fig.  3  by  a  split 
peak  (two  adjacent  peaks).  The  split  peak  which  is  observed  in 
Fig.  3  rear  the  signal  angle  of  +3  degrees  is  apparently  due  to 
interference  of  the  two  noise  peaks  which  are  preser1,  at  angles 
near  the  signal.  Since  the  complex  noise  peaks  are  randomly  located, 
split  peaks  caused  by  such  noise  interference  may  be  eliminated 
with  some  form  of  averaging  within  the  MESA  algorithm.  The  same 
set  of  complex  noise  data  (IR=5)  is  used  in  the  evaluation  of  several 
averaging  techniques . 

V.  AVERAGING  TECHNIQUES 


A.  Averaged  Filter  Weights 

It  is  not  necessary  to  average  a  complete  set  of  filter 
weights.  ,  where  (l^nfN),  since  all  filter  weights  for 

(2in£N)  are  a  function  of  the  last  filter  weight  as  given  by 

eqn.  (3).  The  last  filter  weight,  which  is  given  by  eqn.  (4)  is 


computed  and  averaged  over  L  data  sets. 


N+i 

YN+1 


N+l 

tN+1 


(kAt) 


In  all,  a  total  of  (L-M)  data  points  are  utilized  in  computing  such 
an  averaged  MESA  antenna  pattern. 

An  average  MESA  antenna  pattern  utilizing  an  averaged 
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ignal,  IR- 


last;  filter  weignt  is  shown  in  Fig.  4  for  one  siganl  incident  at 
+5  degrees  and  a  SNR  of  10  dB.  Peak  splitting  is  eliminated  in 
the  example  of  Fig.  \  with  an  average  of  only  two  filter  weights 
(L=2)  and  utilizing  a  total  of  only  16  data  points.  The  16  data 
points  include  the  same  complex  noise  data  set  used  in  the  example 
of  Fig.  2  (IR=5) .  The  signal  is  very  prominent  in  the  averaged 

MESA  pattern  of  Fig.  4,  although  the  signal  peak  is  slightly  displaced 
at  an  angle  of  +4.5  degrees  denoting  an  inaccuracy  of  0.5  degrees. 
Further  averaging  beyond  1=10  provides  little  or  no  improvement. 

The  averaged  MESA  antenna  pattern  for  L=10,  which  is  shown  in  Fig.  5, 
has  nearly  white  noise.  The  noise  peaks  are  very  subdued  and  are 
almost  eliminated,  consequently  very  little  improvement  is  possible. 
But  the  signal  peak,  which  is  considerably  sharpened,  remains  at 
+4.5  degrees  with  an  inaccuracy  of  0.5  degrees. 

Resolution  capability  is  demonstrated  by  the  averaged 
MESA  antenna  pattern  of  Fig.  6,  where  two  signals  with  a  SNR  of 
13  dB  each  signal,  each  element,  are  just  resolved.  Best  resolution, 
which  is  depicted  in  Fig.  6,  is  achieved  for  N=7  and  L=20.  The  two 
signals  as  detected  in  Fig.  6  are  located  closer  together  at  angles 
of  0.5  and  4.5  degrees.  The  technique  of  averaging  filter  weights 
is  a  simple  and  fast  stablization  technique  which  results  in  good 
resolution  and  detection  capability  for  a  relatively  small  number  of 
repetitive  calculations  utilizing  160  (L-M)  data  points. 


B .  Averaged  Prediction  Errors 

An  averaged  .’’ESA  antenna  pattern  may  also  be  computed 
by  averaging  the  forward  and  backward  prediction  errors  as  defined 
by  eqns .  (5)  and  (6).  Prediction  errors  for  filter  sizes  1  -  N 

are  all  calculated  in  the  Burg  technique,  however  best  results 
are  achieved  by  averaging  only  the  set  of  prediction  errors  for 
the  specified  filter  size  N  as  follows: 


for  (NiniM) 


* 
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averaged  filter  weights 


I  ># 


BN  =  i  f  [bN  ,kAt)l  for  ( l£n-M-N ) 

n  B  |  n 

k~l  L  k 

Both  prediction  error  averaging  and  filter  weight  averaging  are 
comparatively  simple  and  fast  averaging  techniques,  but  the  two 
averaged  antenna  patterns  are  quite  dissimilar. 

An  averaged  MESA  antenna  pattern  computed  by  averaging 
the  prediction  errors  is  shown  in  Fig.  7  where  the  split  peak 
(obtained  for  L=l,  Fig.  3)  is  eliminated  by  averaging  Wxth  only 
one  additional  data  set  (L=2) .  The  result  of  further  averaging 
is  shown  for  L=10  in  Fig.  8.  The  filter  size  (N=5)  and  the 

complex  data  set  for  IR=5  are  the  same  as  used  in  all  previous 

examples  of  averaged  MESA  patterns.  It  is  evident  by  observation 
of  Figs.  7  and  8  that  averaging  of  prediction  errors  does  not 
whiten  the  noise  and  does  not  enhance  the  SNR,  but  peak  splitting 

is  eliminated.  In  Fig.  8  the  signal  peak  is  located  at  +4  degrees 

for  an  error  of  one  degree.  Further  averaging  beyond  L=10  does 
not  improve  the  antenna  pattern  for  one  signal  and  an  8  element 
array . 

Prediction  errors  are  averaged  in  Fig.  9  for  two  signals  incident 
at  0  and  +6  degrees.  The  two  signals  are  well  resolved  and  the 
SNR  is  improved  for  L=30.  One  signal  is  accurately  located  at 
0  degrees,  while  the  second  signal  which  is  located  at  +4.5  degrees 
is  in  error  by  1.5  degrees.  The  SNR  is  significantly  improved, 
more  so  than  for  the  single  signal  of  Figs.  7  and  8. 

C .  Averaged  Covariance  Matrix 

While  not  so  obvious,  the  equations  of  the  Burg  technique 

do  contain  elements  of  the  covariance  matrix.  These 

elements  may  be  averaged  and  incorporated  into  the  Burg  technique 

equations  without  altering  the  utilization  or  the  characteristics 

N+ 1 

ol  the  Burg  technique.  The  only  independent  filter  weight 

which  is  defined  by  eqn.  (4)  of  the  Burg  technique,  is  a  function 

N  M 

of  the  forward  and  backward  prediction  errors,  F  and  BT  . 

N  + 1  a. 
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prediction  errors 
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ig.  8  -  One  signal,  averaged  prediction  errors,  t.=l 0 


element 


Products  of  the  prediction  errors  may  be  considered  as  functions 
of  the  covariance  matrix  elements  by  considering  the  original 
definition  as  follows: 


Forward  prediction  error 


FN  0  ?  YN  x 
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Backward  prediction  error 
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where  y^  =  1,0 


The  last  filter  weight  y^J  of  eqn.  (4)  may  be  expressed  as  the 
ratio  of  two  functions,  TOP  and  BOTH,  as  follows: 

N+l 

y,,,7  =  -2  TOP/BOTM 

+  x 
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(10) 


Insertion  of  the  prediction  error  definitions  eqns . 

(7)  and  (8)  into  eqns.  (9)  and  (10)  and  re-ordering  the  summations 
yields  the  desired  functional  form  as  follows: 


TOP  = 


M-N 


1  =  1 


N  N  N  /  N  * 

£  vlN  /  r  (Y  ) 

,,  ,  M  N+I-m+1  “  \  '  n' 

M-- 1  n=  1  v 


Wl 


19 


i 

1 


where 


N  N 

N  N 

E  I  y  Y  r  (N-M-n+2) 
„  ,  ,  m  n 

m= 1  n= 1 


where 


M-N 

r  (N-m-n+2)  =  E  X 


I+N-m+1  I+n-1 


M-N  I  N 


B0™  =  *  E.  Ym  XI+N-m+l  E_  <*“>  XN+ I-n+1 

1=1  I  iti=l  n~l 


N  N  /  i 

N  *  /  N  * 

+  £,  (Ym^  XN+m-l  S  ((Yn5  XI+n-l 

m=l  n=l  V 


NN  NN 

BOTM  -  T,  Z  y”  (y")  r(n-m)  +  l  Z  (y„)  Y„  r(m-n) 


m=l  n=l 


m=  1  n= 1 


M-N 

r(n-m)  =  I  X 


N+I-m+1  N+I  n+1 


M-N 

r(m-n)  =  I  X 


I+m-1  I+n-1 


It  is  apparent  that  the  autocorrelation  coefficients  defined  by 
eqns .  (11) ,  (12)  and  (13)  may  be  averaged  as  follows: 


M-N  I, 

r  (N-m-n+2 )  =  £  -  T.  X(k£t)  X(ki.t) 

I  =  1  JJ  k=l  N+l-m+1  i  +  n-1 
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r  (m-n) 


(15) 


M-N  L  * 

=  Z  ~  T.  X(kAt)  X(kAt) 

1=1  b  k=l  N+I-m+1  N+I-n+1 


M-N  .  1.  * 

r(ir.-n)  ~  Z  ~  Z  X  (kAt)  X  (kAt) 
1=1  L  k=l  I+m-1  I+n-1 


(16) 


The  autocorrelation  coefficients  (eqns.  (11),  (12)  and  (13))  are 

elements  of  the  covariance  matrix  for  the  data  set  X^,  X^,  .  .  .  X^, 

The  averaged  covariance  matrix  elements  given  by  eqns  (14),  (15)  and 

N+l 

(16)  may  be  utilized  to  compute  the  last  filter  weight  and  an 

averaged  MESA  antenna  pattern. 

The  results  due  to  averaging  of  the  covariance  matrix 
are  observed  in  Fig.  10  where  the  split  peak  shown  in  Fig.  3  is 
eliminated  with  only  one  additional  data  set  (L=2).  Also  the  noise 
is  considerably  whiter  which  greatly  improves  the  SNR.  Further 
averaging  provides  little  additional  improvement  as  noted  by  Fig.  11, 
where  for  L=10  the  noise  appears  slightly  whiter  and  the  SNR  is 
slightly  enhanced  over  the  results  shown  in  Fig.  10. 

Two  signals  located  at  0  and  +6  degrees  are  well  resolved 
in  Fig.  12  with  covariance  matrix  averaging,  for  L=10  and  an  input 
SNR  of  13  dB  each  signal.  The  background  noise  is  substantially 
reduced  although  it  is  not  as  white  as  the  noise  that  appears  in 
Fig.  11.  Averaging  of  the  covariance  matrix  elements  is,  as  demon¬ 
strated,  an  excellent  averaging  and  stabilization  technique. 

D.  Averaged  MESA  Snapshots 

Individual  MESA  antenna  patterns,  which  are  referred  to 
as  "snapshots"  (e.g.  Fig.  la,  lb),  may  be.  averaged  in  order  to 
obtain  a  stable  antenna  pattern.  This  method  has  been  utilized  (14) 
successfully  using  optimal  filter  sizes.  However,  if  the  optimal 
filter  size  cannot  be  determined,  a  fixed  filter  size  may  be  selected 
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Fig.  10  -  One  signal,  averaged  covariance  matrix. 
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for  computing  all  the  individual  MESA  snapshot  antenna  patterns. 

The  best  resolution  is  always  achieved  using  the  largest  possible 
filter  size  (N=M-1).  However,  the  larger  filter  sizes  are  also 
the  most  unstable.  Consequently,  while  an  average  of  many  MESA 
snapshots  improves  stability,  the  resultant  stable  pattern  may  not 
be  a  very  accurate  antenna  pattern. 

Averaging  results  for  two  filter  sizes  (N=5  and  N=7) 
are  demonstrated  in  the  three  following  examples.  A  simple  average 
of  two  MESA  snapshots  (L=2)  is  shown  in  Fig.  13  for  one  signal 
incident  at  5  degrees  with  a  SNR  of  10  dB.  The  split  peak  that 
occurs  for  L=l  (Fig.  3)  is  still  present  in  the  average  of  two 
snapshots.  However,  there  is  improvement  in  the  SNR.  There  are 
cf  course  twice  the  number  of  peaks  (10)  as  expected  for  two  MESA 
snapshots  having  five  filter  weights  (N=5)  each.  The  consequence 
of  further  averaging  is  demonstrated  in  Fig.  14  where  30  antenna 
patterns  (computed  for  a  signal  incident  at  5  degrees  with  a  SNR 
of  10  dB)  are  averaged.  There  is  further  improvement  in  the  signal 
peak  definition  and  accuracy;  the  signal  is  located  at  +4  degrees, 
for  an  error  of  one  degree.  The  SNR  is  improved  substantially . 

In  Fig.  15  the  resultant  "'nrago  of  MESA  antenna  patterns  for  two 

signals  incident  at  0  and  T6  degrees  is  disappointing  as  the  two 
signal  peaks  are  not  very  well  defined.  Instead  there  are  four 
strong  peaks,  two  of  which  are  in  error.  There  is  of  course 
improvement  in  the  SNR,  but  the  resolution  characteristics  are  very 
poor.  In  the  example  of  Fig.  15,  averaged  MESA  snapshots  ail  have 
the  maximum  number  (7)  of  filter  weights.  Better  results  have 
been  achieved  (14)  using  an  optimal  filter  size,  however  the 
optimal  filter  size  can  only  be  computed  if  the  incident  signal 
angles  are  known. 

VI.  ADAPTIVE  FILTER  WEIGHTS 

While  the  MESA  technique  is  inherently  adaptive,  other 
adaptive  methods  which  have  been  demonstrated  (17,  18)  have  simple 
procedures  for  updating  the  filter  weights.  One  such  procedure  (12) 
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increases  the  filter  weights  in  proportion  t.o  the  correlation  of 
the  prediction  error  with  the  data  set.  For  example,  a  filter 
weight  W(k)  computed  at  time  kAt  may  be  updated  at  a  later  time 
(k+1)  At  to  give  W(k+1)  as  follows: 


W (k+1) 


W  (k) 


+ 


M 

V 

n=l 


en(k)X*(k) 


(17) 


where  the  correlation  is  taken  over  all  computed  prediction  errors 

e  (k) .  The  prediction  errors  are  defined  as  follows: 
n 


en(k)  -  Xn(k>  -  £n(k)  (18) 

for  data  samples  Xn(k)  and  predicted  values  ”x n  ( k ).  The.  proportionality 
constant  is  denoted  by  the  convergence  parameter  y . 

In  order  to  incoporate  this  procedure  into  uhe  MESA 

technique,  a  first  set  of  filter  weights,  is  computed  in  the  usual 

manner  as  defined  by  eqns.  (3)  and  (4).  Subsequent  filter  weights 
N 

r  (k+1)  mc.y  then  be  computed  according  to  eqn.  (17)  as  follows: 


,N 


n 


(k+1) 


N 

r”  (k)  +  y 


M-N  * 

1  e”(k)X  (k) 
n=l  11  n 


(19) 


where  the  prediction  error  is  actually  the  sum  of  the  forward  and 
backward  predic  ion  errors  over  all  poss'blo  errors  as  follows: 

4{k)  =  Fn+W(k)  f  Bn(k)  (20) 

However,  with  use  of  the  Burg  technique,  onlv  the  last  filter  weight 

„N 

,,j  ,  need  be  computed  with  eqn.  (19),  since  ail  other  filter  w-iqhts 
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(n<N)  are  dependent  upon  according  to  eqn.  (3). 

As  the  prediction  error  is  whitened,  the  additive  adaptive 
component  is  reduced,  since  the  correlation  of  a  whiter  prediction 
error  with  the  data  set  is  smaller.  Consequently,  the  adaptive 
filter  weights  may  converge  to  become  a  whitening  filter. 

The  result  of  updating  MESA  filter  weights  is  illustrated 
in  Fig.  16  where  an  "adapted"  MESA  antenna  pattern  is  hown  for  one 
signal  incident  at  +5  degrees,  SNR=10  dB,  and  L,=  2.  The  "adapted" 
pattern  for  L=2  is  of  course  quite  similar  to  the  computed  MESA 
snapshot  (L=l)  shown  in  Fia.  3,  since  the  filter  weights  have  been 
modified  only  once.  The  split  peak  is  still  present  for  L=2  in 
Fig.  16,  but  with  further  adaption  the  split  peak  is  eliminated  as 
shown  in  Fig.  17  for  L=10.  However,  the  SNR  is  not  improved,  although 
three  noise  peaks  have  been  reduced.  As  observed  in  Fog.  17  the 
noise  has  remained  peaked  even  after  ten  adaptions.  Obviously  the 
adaption  method  does  not  tend  to  whiten  the  noise,  and  consequently 
the  results  are  most  disappointing. 

In  another  application  of  the  adaptive  method,  two 
signals  incident  at  0  and  +6  degrees  are  resolved  in  con  adaption 
(L=10)  as  shown  in  Fig.  18,  where  the  SNP  is  13  dB  for  each  signal. 

One  signal  is  located  accurately  at  +6  degrees  while  the  other 
signal  is  located  at  +1.?  degrees  with  an  error  of  1.5  degrees. 

The  SNR  is  improved  with  respect  to  the  original  (L=l)  MESA 
snapshot  of  Fig.  3.  However,  the  noise  in  not  whitened  in  the 
adaptive  process  as  had  been  anticipated. 

The  results  from  using  adaptive  filter  weights  with 
MESA  are  very  disappointing,  and  in  addition  the  value  of  a  con¬ 
vergence  parameter  must  be  specified.  In  the  three  preceding 
examples  of  the  adaptive  technique  the  value  of  the  convergence 
parameter  y  is  quite  critica.  If  v  is  too  small,  there  is  little 
improvement  in  the  computed  antenna  pattern,  and  if  y  is  too  large, 
there  may  be  considerable  distortion  of  the  signal  peak. 
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adaptive 


signal,  adaptive  filter  weights  L=10 


VII.  RESOLUTION 


Of  the  five  stablization  methods  examined,  only  two 
appear  to  preserve  and  enhance  the  desirable  MESA  characteristic 
of  whitening  the  noise.  The  resolution  capability  of  both  whitening, 
stablization  methods  is  demonstrated  in  one  example  of  two  signals 
separate  l  by  four  degrees,  at  -2  and  +2  degrees,  with  equal  SNR 
values  of  20  dB  each  antenna  element. 

A.  st.ablized  antenna  pattern  computed  by  averaging  the 
covariance  matrix  is  shown  in  Fig.  19,  where  the  two  signals  are 
not  resolved,  but  only  one  signal  peak  is  detected  at  +0.5  degrees. 
However  both  signals  are  identified  in  Fig.  20,  where  the  antenna 
pattern  is  computed  by  averaging  the  filter  weights  using  a  short 
average  of  L=5.  The  same  short  average  is  used  in  computing  both 
antenna  patterns  in  Figs.  19  and  20  since  further  averag.imr,  which 
does  appear  to  improve  signal  detection,  only  serves  to  reduce  the 
resolution  capability.  However,  some  averaging  if  necessary  in  order 
to  obtain  a  stable,  reliable  antenna  pattern.  The  noise  is  whiter 
in  Fig.  19,  but  the  resolution  Is  best  in  Fig.  20.  While  the  two 
character istics  appear  to  be  somewhat  incompatible,  it  is  necessary 
that  good  resolution  be  achieved  with  stable,  reliable  antenna 
patterns . 
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Two  unresolved  signals,  averaged  covariance  matrix 
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VIII.  CONCLUSIONS 


Split  signal  peaks,  which  are  a  common  occurence  in 
MESA  snapshot  patterns,  are  shown  to  be  a  consequence  of  noise 
interference.  Such  noise  interference  is  virtually  eliminated 
with  use  of  the  proposed  stablization  methods. 

Of  the  five  stablization  techniques  examined,  two  have 
excellent  characteristics,  one  other  is  only  somewhat  satisfactory, 
and  two  were  very  disappointing.  An  average  of  filter  weights 
and  an  average  of  the  covariance  matrix  are  both  very  useful 
stablization  methods.  Both  methods  serve  to  whiten  the  noise 
and  greatly  improve  the  SNR.  In  addition,  split  signal  peaks 
were  not  observed  with  use  of  either  aveiaged  filter  weights 
or  averaged  covariance  matrix  elements.  Further  testing  of 
these  two  averaging  methods  is  clearly  justified.  Hopefully, 
resolution  and  SNR  properties  of  these  two  excellent  averaging 
and  stablization  techniques:  will  be  specifically  determined  in 
future  research  efforts. 

It  is  doubtful  that  any  of  the  other  three  examined 
stablization  methods  are  worthy  of  further  consideration. 

Neither  the  averaged  prediction  errors  nor  the  adaptive  filter 
weights  served  to  whiten  the  noise,  and  the  averaged  MF.SA  patterns 
proved  to  be  most  unstable. 

In  the  one  example  of  two  signals  separated  four 
degrees  the  MESA  antenna  pattern  computed  by  averaging  filter 
weights  provided  the  best  signal  resolution.  However,  many  such 
examples  need  to  be  accumulated  in  order  to  determine  the  resolution 
characteristics  of  the  recommended  MESA  stablization  methods. 
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BROADBAND  SIGNAL  DETECTION  WITH  TIIE  SCOT  AND  SCOT-PET 

1,0  INTRODUCTION 

1 . 1  Background 

The  crosscorreiatiori  technique  is  useful  for 
detecting  broadband  radar  pulses.  However,  it  is  well 
known  that  the  crosscorrelation  technique  is  adversely 
affected  by  narrowband  interference  or  narrowband  signals 
transmitted  by  radar  surveillance  systems.  Consequently, 
the  crosscorrelation  peak  (due  to  the  broadband  signal) 
may  be  obscured  by  a  sinusoidal  modulation  of  the  cross¬ 
correlation  function  caused  by  the  presence  of  strong 
narrowband  components. 

Other  si.grial  processing  techniques  may  not  be 
so  adversely  affected  by  the  presence  of  narrowband  signal 
components.  Carter,  et  al.  (1)  have  provided  one  snapshot 
example,  where  a  broadband  signal  in  the  presence  of  three 
nt rrowband  signals  is  not  detected  by  the  crosscorrelation 
function,  but  is  readily  detected  with  the.  smoothed  coher¬ 
ence  transform  (SCOT) . 

The  SCOT  is  the  Fourier  transform  of  a  cross¬ 
power  spectral  function,  which  has  ?  frequency  dependent 
normalization.  The  normalization  serves  to  whiten  the 
crosspower  spectra  and  thereby  minimize  the  effect  of 
narrowband  signal  interference.  While  the  SCOT  may  pro¬ 
vide  improved  broadband  signal  detection,  it  may  be  possible 
to  improve  the  SCOT,  as  defined  by  Carter,,  et  al.,  by  esti¬ 
mating  the  crosspower  function  with  use  of  the  prediciton 
error  transform  (PET).  King  (2)  has  demonstrated  that  the 


autocorrelation  function  of  broadband  signals  may  be  more 
accurately  evaluated  using  the  PET. 


1 . 2  Comparative  Evaluation 

In  order  to  compare  and  evaluate  the  SCOT,  SCOT- 
PET,  and  the  cross-correlation  function,  a  set  of  receiver 
operating  characteristic  (ROC)  curves  are  computed  using 
data  simulated  for  two  time  dependent  function  x(t)  .and  y(t).  For 
each  function  the  simulated  data  contains  a  coherent,  white, 
broadband  signal,  incoherent,  white  oroadband  noise,  and 
four  coherent,  narrowband  non-harmonic  s  i.gnals . 

2.0  THEORY 


2  .  .1  Cross  correlation 

The  crosscorrelation  function  may  be  defined  in 
the  time  domain  by  a  convolution  integral,  or  in  the  fre¬ 
quency  domain  by  the  Fourier  transform.  Since  it  is  con¬ 
venient  in  this  investigation  to  simulate  sensor  data  in 
the  frequency  domain,  the  crosscorrelation  function  is 
defined  as  follows: 


^xy(T)  =  [ 


Gxy(f)  ei2*1  df 


:d 


where  G  (£)  is  the  crosspower  spectral  density  function, 
xy 

The  crosscorrelation  may  be  normalized  by  the  following 
equation : 
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where  <t>  (t)  and  <J>  (t)  are  the  autocorrelation  functions  for 

*  y 

the  respective  signals  x(t)  and  y(t)  and  G^tf)  and  G^(f)  are 
the  auto  power  spectra  evaluated  for  the  respective  signals 
x(t)  and  y(t).  The  normalised  crosscorrelation  as  defined 
by  eqn.  (2)  is  evaluated  for  comparison  with  the  SCOT,  which 
is  also  a  normalized  function. 


2 . 2  Smoothed  Coherence  Transform  (SCOT) 

The  SCOT  is  defined  by  Carter  et  al.  (1)  as 

follows : 


C(x)  =  /  W(f  )y  (f)  el27Tfldf  (5) 

—  oo 

where  W(f)  is  a  weighting  function  of  choice,  and  Y(f) 
is  the  crosspower  spectral  density  function  with  a  frequency 
dependent  normalization.  The  spectral  function  is  defined 
as  follows: 


xy  (f) 

Y(f)  = - ===== 

/G  (f)G  (f) 

x  y 


(6) 


where  G  (f)  ,  G  (f)  and  G  (f)  are:  the  crosspower  and  auto- 
xy  x  y 

power  spectral  density  functions  evaluated  for  the  respective 


signals  x(t)  and  y(t).  For  the  applications  considered  in 
this  investigation  all  power  spectral  density  functions  are 
defined  in  the  frequency  domain. 

2 . 3  The  SCOT-PET  Function 

The  autocorrelation  function  may  be  represented 
by  the  prediction  error  transform  (PET)  as  demonstrated  by 
King  (2)  and  in  a  similar  manner  the  crosscorrelation  and 
SCOT  functions  may  also  be  represented  by  PET.  Actually, 

PET  is  simply  an  inverse  Fourier  transform,  which  has  a. 
derivation  based  upon  a  prediction  error  function.  When 
applied  tc  SCOT,  the  cross -spectral  function  y(f)  is  esti¬ 
mated  or  predicted  with  a  discrete  convolution  filter,  and 
a  prediction  error  is  defined  in  the  frequency  domain. 

The  maximum  entropy  power  spectra  has  a  similar  PET  repre¬ 
sentation  based  upon  a  prediction  error  defined  in  the 
time  domain. 

The  PET  may  be  a  useful  representation  of  a  function 
that  is  derivable1  from  a  Fourier  transformation.  The  PET 
representation  is  useful  only  if  the  transformed  function 
(such  as  y(f))is  known  over  a  limited  region  in  the  time 
or  frequency  domains.  in  order  to  derive  the  PET  the 
transformed  function  must  also  be  a  predictable  (non-random) 
function  or  have  predictable  components. 

The  cross-power  spectral  function  y(t)  has  a 
predictable,  periodic  modulation  given  by 

i  2irfT 
e 

where  T  is  the  time  delay  between  two  broadband  signals 
X(t)  and  X(t+T)  which  have  a  common,  coherent,  broadband 
component.  If  the  periodic  modulation  is  well  defined 
in  a  low  noise  environment  over  a  sufficient  spectral 
interval,  then  the  Fourier  transform  of  the  modulating 
component,  exp(i2irfT)  ,  is  a  well  defined,  detectable  SIN'’ 
function  centered  about  T  in  the  time  domain.  However, 
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if  noise  or  other  interfering  components  are  present  in  the 
cross  power  spectral  function  y(f)  ,  then  the  SCOT  function, 
which  is  defined  by  the  Fourier  Transform,  may  not  provide 
a  well  defined  SINC  function  representative  of  the  coherent, 
broadband  signal  component.  When  interference  or  noisy 
components  are  present,  or  when  the  cross-power  spectral 
function  y(f)  is  known  only  over  a  limited  bandwidth,  then 
the  periodic  modulation  component  of  y (f)  may  be  better 
defined  with  use  of  the  prediction  filter  a  as  follows: 


T.  c 

n=l 


f-n 


(7) 


where  y^  is  the  discrete  representation  (N  components)  of 
the  function  y(f).  A  prediction  error  may  be  defined  by 
the  expression 


=  Y, 


!3) 


N 

and  a  new  filter  b^,  the  prediction  error  filter,  may  be 
j.ntroducted  as  follows: 


ef 


N 

T.  b 
n=0 


N 

n 


f-n 


(9) 


If  it  is  recalled  that  the  SCOT,  C(i)  ,  is  defined  as  the 
Fourier  transform  of  y(f)  (  then  the  Fourier  transform  of 
egn.  (9)  results  in  the  following  expression: 


E_. 

IN 


C  (t) 


N 

1  b 
n=0 


N  -i2irtn  ( Af ) 
e 
n 


(10) 


where  B.,  is  the  Fourier  transform  of  ec.  The  SCOT  function 
N  r 

is  defined  by  solving  eqn.  (10)  as  follows: 
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C  (  T  } 


(11) 
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N 

E 

n=0 


-i2irtn  ( Af ) 
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The  SCOT,  as  defined  by  eqn.  (11),  is  the  inverse  of  the 

Fourier  transform  of  the  prediction  error  filter  b£j.  The 

inverse  Fourier  transform  solution  is  called  the  prediction 

error  transform,  because  its  derivation  is  a  result  of  the 

definition  and  application  of  the  prediction  error  filter. 

N 

The  prediction  error  filter  coefficients,  b  , 
r  n 

may  be  evaluated  by  squaring  and  minimising  the  error 

given  by  eqn.  (9).  The  solutions  for  the  coefficients 
N 

b  and  the  constant  E„  have  been  formulated  by  King  (3) 
n  N 

with  use  of  the  Burg  technique.  The  solution  agrees  with 

that  obtained  by  Burg  (4) ,  who  found  the  unknown  coefficients 

N 

by  maximizing  the  entropy.  solutions  for  and  b^  are  a 
set  of  iterative  equations  listed  as  follows: 


E.  =  r 
1  o 


(12a) 


(r  ■-  zero  delay  autocorrelation  coefficient 
o 

of  data  set  consisting  of  M  data  samples) 
2 


1,-1  -  eh  [ 


1  *  <bsti>r 


b"  =  1.0 


(12b) 

(13a) 
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,  N+l 
b..._  = 
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,„N  *  N 
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E 
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(13b) 


,  0N ,  “  ,  ,  N  S' 

<V  +  a  j  +N  j 


,N+ 1  ,  N  A  .N+l  ,,N  ' 

bn  =  bn  +  bN+l  bN-n+2 
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The  forward  prediction  errors 
'  M 

tion  errors  S  .  are  defined  as 
J 

power  spectral  data  set  as 


,  and  the  backward  prcdic- 
a  function  of  the  known,  cross- 
f ollows : 


(14c.) 


„  N+l 
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3-N 


8N+1 
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/kN+1. 
bN+ 1 


N 
a.  . 
j+M 


(14b) 

(15a) 

(15b) 


3.0  DATA  SIMULATION 


Since  the  crosscorrelation,  SCOT,  and  SCOT-PET 
functions  must  necessarily  be  evaluated  in  succession  and 
repeatedly  in  order  to  plot  a  set  of  ROC  curves,  it  is 
wise  to  simulate  data  in  the  frequency  domain  if  possible, 
so  as  to  prevent  repeated  transformations  of  data  sets 
otherwise  simulated  in  the  time  domain.  Both  the  cross¬ 
correlation  and  SCOT  functions  require  a  representation 

of  the  c'i  osspower  spectral  function  0  (f).  The  spectral 

xy 

function  C  (f)  may  be  expressed  very  simply  by  assuming 
that  the  broadband  and  narrowband  signal  components  are 
steady  state  signals  such  that 


Gxy(f) 


Vf) 


i2ir£T 


where  G„ (f)  is  the  an to power  spectra  of  the  coherent  br<  ad- 
bard  component  and  G„  is  the  aulowwei  of  the  coherent  narrow- 

*•  ;n 

band  component  and  the  time  delay  between  the  funcions  x(t) 
and  y ( L)  is  T  seconds.  Since  it  is  assumed  that  the  signals 
• :  ( t )  and  Y ( t )  both  contain  incoherent,  broadband  components. 
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the  respective  autopower  spectral  functions  G^(f)  and 
Gy(f)  are  given  by: 


Gx(f)  =  GI(t)  +  GB(f)  +  GN 


(17) 


Gy(f)  =  Gx(f)  +  GB(f)  +  Gn 


where  the  incoherent  broadband  components  G^.(f)  is  assumed 
to  be  of  equal  power  for  both  x(t)  and  y(t) . 

The  signal  component  spectra  may  take  any  form 
over  a  given  bandwidth,  but  in  order  to  establish  a  standard 
for  future  ROC  curve  comparisons,  all  broadband  spectra 
(incoherent  and  coherent)  are  defined  to  be  uniform  (white; 
spectra  over  an  arbitrary  bandwidth.  By  assuming  white 
broadband  spectra,  the  whitening  effect  of  frequency  depen¬ 
dent  normalization  (used  in  the  SCOT  function)  is  minimized. 
Consequently,  the  improvement  in  signal  detection,  provided 
by  the  SCOT,  is  a  minimal  improvement  due  only  to  the 
normal:. zation  of  the  narrowband  components. 

The  broadband  spectral  components  are  assumed 
to  have  signal  amplitudes  with  Gaussian  distributions  such 
that, 


Gx(f) 

=  „2 
'  I 

In  ( 1/Rj ) 

(18) 

Vf) 

2 

=  °B 

ln(l/RB) 

(19) 

whe  :e  1^  and  RQ  are  uniformly  distributed  random  numbers. 
Values  are  assigned  by  arbitrari  ly  let-Hng 

oT  =  1.0 

and  by  defining  a  broadband  power  SNR  as  follows: 

BBSNR  -  10  bog (o^/ol) 

J. 

a  ( 

BBSNR  =  10  Log(O^) 


4  •****-?"•  . 


Similarly,  a  narrowband  SNR  is  defined  with 
respect  to  unit  power  so  that, 

NBSNR  =  10.  Log(G  )  . 


(21) 


It  remains  only  to  specify  values  of  the  cross¬ 
correlation  normalization  functions  <t>x(o)  and  ^(o)  as  follows: 

<f>x(o)  =  /  [a1(f)  +  gb  (f )]  df  +  gn 

(22) 


4>„(o)  =  2  G_  +  G 


[< 


B 


+  G. 


N, 


where  the  spectral  components  are  arbitrarily  defined  to 
exist  over  a  spectral  band  from  -1  to  +1.  Since  identical 
component  power  levels  are  assumed  in  both  signals  x(t.)  and 
Y  (t)  , 


<J>  (o)  =  4>  (o)  (2  3) 

y  x 

The  crosspower  spectral  functions  (f)  ,  <j>  (f) 

and  the  autopower  functions  G  (f)  ,  G  (f),  6  (o',  <t>  (o)  are 

x  y  x  / 

all  specified  by  the  preceding  equations  for  given  values 
of  BBENR  and  NBSNR.  The  crosscorrelation,  SCOT,  and  SCOT- 
PET  are  Fourier  or  Prediction  Error  transforms  of  the 
defined  crosspower  and  autopower  spectral  functions  as 
defined  by  eqns .  (2),  (5),  (6)  and  (11).  A  set  of  ROC 

curves  may  be  constructed  by  re-evaluutiny  the  cross¬ 
correlation,  SCOT,  and  SCOT-PET  functions  repeatedly,  and 
counting  the  false-alarm  peaks  and  signal  peaks  above 
specific  threshold  values. 

4.0  EXAMPLES  OF  CROSSPOWER  FUNCTIONS 


The  complex,  crosspower  spectral  function  y  (f) 
is  computed  and  plotted  in  Figs.  1,  2  and  3.  The  spectral 
function  consists  of  white,  broadband,  coherent  and  incoherent 
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components  and  four  equal  power  narrowband  components.  The 
magnitude  function  f y ( f ) j  is  observed  in  Fig.  1  for  128 
computed  points,  and  for  a  broadband  SNR  =  1. 55  dB ,  and  a 
narrowband  SNR  =  20  dB .  The  real  parr  of  y  (f)  is  shown 
for  128  computed  points  and  the  same  SNR  values  in  Fig.  2. 
The  real  part  of  y  (f)  is  also  shown  in  Fig.  3  for  123 
computed  points,  but  with  the  BBSNR  =  -20  dB  and  the 
NBSNR  =  10  dB. 


The  low  frequency  modulation  observed  in  Fig.  2 

is  not.  real,  but  is  rather  due  to  a  sampling  rate  which  is 

inadequate  for  display  purposes,  although  it  is  adequate 

for  computational  purposes.  The  higher,  frequency  component 

i2irfT 

e 

is  well  defined  in  Fig.  2,  but  not  so  visible  in  Fig.  3 
due  to  the  large  difference  in  BBSNR.  Of  course  when  the 
B8SNR  is  reduced  as  in  Fig.  3,  the  tour  narrowband  compo¬ 
nents  are  more  visible. 


5.0  EXAMPLES  OF 

CROSSCORRELATION,  SCOT,  AND  SCOT-PET  TIME  FUNCTIONS 

The  normalized  crosscorrelator  time  function 
(in  decibel  units)  is  shown  in  Fig.  4  for  a  time  delay  of 
-0,3  sec.,  BBSNR  =  -20  dB,  and  four  narrowband  components 
of  10  dB  each.  The  broadband  signal  peak,  is  visible  at  the 
time  delay  of  -0.3  sec.,  but  only  slightly  above  the 
largest  clutter  peaks.  The  SCOT  time  function  (in  decibel 
units)  computed  for  the  same  SNR  values  (as  used  in  Fig.  4) 
has  a  more  prominent  signal  peak  as  observed  in  Fig.  5. 

The  signal  peak  height  is  the  only  notifiable  difference 
between  the  crosscorrelator  and  SCOT  time  functions.  The 
clutter  patterns  for  each  time  function  are  identical, 
because  the  broadband  signal  and  noise  spectra  are  uniform 
(white)  spectral  distributions  The  clutter  patterns  would 
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be  very  different  if  the  broadband  spectra  were  not  uniformly 
(non-white)  distributed,  because  the  SCOT  has  a  frequency 
dependent  normalization.  Consequently,  the  detection  charac¬ 
teristics  of  the  SCOT  and  SCOT-PET  functions  are  expected 
to  improve  significantly  for  non-white  broadband  spectral 
functions . 

The  SCOT-PET  time  function,  which  is  computed  for 
the  same  SNR  parameters  used  in  Figs.  4  and  5,  is  shown 
plotted  in  decibels  in  Fig.  6.  The  signal  peak  is  even 
more  prominent  at  the  time  delay  of  -0.3  sec.,  and  fewer 
clutter  peaks  are  observed  than  in  the  two  previous  time 
functions.  Based  upon  the  three  detector  time  functions 
shown  in  Figs.  4-6,  the  SCOT-PET  ti.me  function  appears 
to  offer  the  best  signal  detection  capability. 

6.0  RCC  CURVES 

6 . 1  Construction Based  Upon  100  Time  Functions 

One  set  of  time  functions  is  insufficient  evi¬ 
dence  to  form  judgment,  so  a  set  of  receiver  operating 
characteristic  (ROC)  curves  are  constructed  based  upon 
100  sots  of  computed  time  functions.  The  results  are 
shown  in  Figs.  7,  8  and  9.  The  ROC  curves  are  constructed 
for  false  alarm  probabilities  between  .01  and  1.0.  Lower 
values  of  the  false  alarm  probability  would  require  that 
many  more  correlation,  SCOT,  and  SCOT-PET  time  functions 
would  necessarily  be  computed.  While  the  range  of  false 
alarm  probabilities  investigated  are  not  necessarily  of 
practical  value,  they  are  adequate  for  comparing  the  three 
different  detectors.  The  set  of  three  ROC  curves  are  con¬ 
structed  for  a  BBSNF  =  -25  dB  and  a  NBSNR  =  10  dB  (for  each 
of  four  MB  signals) ,  The  crosscorrelation,  SCOT,  and  SCOT- 
PET  ROC  curves  are  shown  in  Figs.  7,  8  and  9  respectively. 

Of  the  three  ROC  curves  shown  the  SCOT  has  the  best  detection 
probabilities.  Of  course  the  lower  false  alarm  probabilities 
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are  usually  of  more  interest  in  underwater  acoustics  appli¬ 
cations,  but  only  for  detection  probabilites  above  0.5. 

The  3C0T-PET  ROC  curve  has  the  lowest  false  alarm  proba¬ 
bility  for  the  detection  probability  of  0.5. 

The  constructed  ROC  curves  are  more  reliable  in 
regions  of  high  false  alarm  probabilities  due  to  the  high 
density  of  data  collected  in  those  regions.  Consider  that 
data  points  with  a  detection  probability  of  1.0  are  based 
upon  100  signal  detections  whereas  a  data  point  representing 
a  detection  prob.ibility  of  0.01  is  based  upon  only  one 
detected  signal. 


6.2  Other  BBSNR  Values 


Two  other  sets  of  ROC  curves  are  constructed  for 
BBSNR  values  of  -40  dB  and  -20  dB  to  illustrate  the  similari¬ 
ties  and  differences  of  the  three  broadband  signal  detectors 
with  changing  BBSNR  values.  Again,  four  narrowband  signals 
having  NBSNR  =  10  dB  each  are  employed  in  constructing  the 
ROC  curves.  In  order  to  minimize  computing  time  and  costs, 
these  and  other  ROC  curves  are  constructed  from  only  25 
sets  cf  time  functions.  However,  these  ROC  curves  are  pre¬ 
sented  only  to  indicate  trends  due  to  parameter  variation. 

In  Fig.  It  the  crosscorrelator,  SCOT  and  SCOT-PET  ROC  curves 
are  shown  on  the  same  graph  for  a  BBSNR  =  -40  dB.  The  cross¬ 
correlator  and  SCOT  ROC  curves  are  nearly  identical  for  all 
false  alarm  probabilities,  whereas  the  SCOT-PET  ROC  curve 
has  shifted  toward  higher  false  alarms.  Detection  perform¬ 
ance  has  deteriorated  significantly  for  the  SCOT  and  even 
more  so  for  the  SCOT-PET.  Ho.rever,  the  detection  perform¬ 
ance  has  changed  only  slightly  for  the  crosscorrelator. 
Apparently,  the  SCOT  approaches  the  crosscorrelator  perform¬ 
ance  and  the  SCOT-PET  performance  deteriorates  toward 
lower  detection  probabilities  with  decreasing  BBSNR  values. 

The  set  of  ROC  curves  constructed  for  BBSNR  =  -20  dB 
is  shown  in  Fig.  11.  The  SCOT-PET  ROC  curve  clearly  indicates 


shown  in  Fig.  11. 


superior  detection  performance  at  lower  false  alarm  proba¬ 
bilities,  although  the  SCOT  is  also  considerably  better  than 
the  crosscorrelator  for  lower  false  alarm  probabilities. 

Again  the  crosscorrelator  ROC  curve  has  changed  only  slightly 
toward  an  improved  detection  performance.  It  is  apparent 
that  for  the  BBSNR  values  investigated,  the  crosscorrelator 
detection  performance  is  inadequate  for  applications  requiring 
a  low  false  alarm  probability.  But  both  the  SCOT  and  SCOT- 
PET  may  have  useful  detection  characteristics  in  regions  of 
low  false  alarm  probabilities  for  EB3NR  values  of  -25  dB  or 
lower.  However,  this  investigation  does  not  examine  detec¬ 
tion  performance  capability  in  regions  of  low  false  alarms. 

6 . 3  Partial ly  Coherent  Broadband  Noise 

All  previous  discussions  and  examples  have  assumed 
that  the  broadband  noise  (occurring  in  the  signals  x(t)  and 
y(t)  is  incoherent.  But  complete  incoherence  is  unlikly  in 
most  radar  clutter  environments.  For  example,  low  level 
partially  coherent,  directional,  broadband  signals  may  be 
present  in  the  atmospheric  clutter  or  in  receiver  channels.  There¬ 
fore,  a  set  or  ROC  curves  are  shown  in  Fig.  12  for  a  3BSMR  =  -25  dB 
and  four  narrowband  components  having  a  NBSNR  -  10  dB  each.  The 
ROC  curves  of  Fig.  12  indicate  that  the  detection  performance 
of  all  three  detectors  has  deteriorated,  but  by  very  different 
amounts.  The  SCOT-PET  has  suffered  the  most  pprfnrmsnce  deteri¬ 
oration,  while  the  crosscorrelator  has  only  suffered  a  slight 
deterioration.  While  the  SCOT  detection  performance  is  dimin¬ 
ished  by  the  noise  partial  coherence,  the  SCOT  detection  perform¬ 
ance  remains  superior  to  that  of  the  crosscorrelator. 

7 . 0  SUMMARY 

The  broadband  signal  detection  performance  of  the 
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crosscorrelator,  SCOT,  and  the  SCOT-PET  have  been  compared 
for  an  environment  containing  four  strong,  interfering, 
narrowband  signal  components,  and  strong,  independent, 
broadband  noise.  Under  such  conditions,  the  crosscorrelator 
e:ihibits  poor,  but  stable  detection  performance,  whereas 
the  SCOT  and  SCOT-PET  appear  to  have  useful  detection  charac¬ 
teristics  in  regions  of  low  false  alarm  probabilities. 
However,  the  SCOT  and  SCOT-PET  detection  performance  is 
inhibited  by  increasing  noise  levels  and  increasing  noise 
coherence.  The  SCOT-PET  has  the  best  detection  performance 
for  low  and  incoherent  noise  conditions,  but  the  SCOT 
detection  performance  becomes  superior  as  noise  levels  and 
noise  coherence  increases. 

Unfortunately,  the  constructed  ROC  curves  do 
not  indicate  the  actual  detection  performance  of  the  SCOT 
and  SCOT-PET  in  regions  of  very  low  false  alarm  probabilities 
Most  detector  applications  require  very  low  false  alarm  proba 
bi litres.  However,  the  constructed  ROC  curves  do  show  that 
the  SCOT  and  SCOT-PET  have  application  in  the  presence  of 
strong,  narrowband  signal  components,  which  very  severely 
degrade  the  detection  performance  of  the  crosscorrelator. 

Both  the  SCOT  and  SCOT-PET  will  have  even  better 
detection  capability  in  the  presence  of  non-white,  b>  oad 
band  noise.  Also,  the  detection  performance  of  the  SCOT- 
PET  may  be  improved  further  with  use  of  larger  filter  sizes 
and  when  fewer  cycles  of  the  modulation  component 

j.wt 

e 

are  present.  These  initial  conclusions  indicate  that  the 
SCOT  and  SCOT-PET  deserve  serious  consideration  as  broad¬ 
band  signal  detectors,  and  that  their  detection  character¬ 
istics  should  be  investigated  further. 
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ROC  CURVES 


ALARM  PROBABtUm 


RESOLUTION  AND  DETECTION  CHARACTERISTICS 
OF  THE  MAXIMUM  ENTROPY  METHOD 


1 . 0  Introduction 

Investigations  of  the  maximum  entropy  method  (MEM) 
nearly  always  provide  an  examination  of  the  technique  and  its 
properties  with  the  use  of  simulated  data.  It  is  thought 
that  this  report  describes  one  of  the  earliest  investigations 
of  MEM  using  actual  radar  signals  that  are  received  with  a 
uniform,  linear,  antenna  array.  Stable  MEM  wavenumber  spectra 
that  are  computed  using  the  actual  data  apparently  have  the 
same  high  resolution  that  has  been  so  often  demonstrated  with 
simulated  data.  Simulated  data  cannot  include  all  the  charac¬ 
teristics  of  actual  radar  data,  which  is  usually  strongly  ani¬ 
sotropic  in  wavenumber  and  time  variant.  Therefore,  tests 
of  the  MEM  that  use  actual  radar  data  are  most  significant  in 
demonstrating  the  usefulness  and  accuracy  of  MEM. 

Since  the  MEM  is  effective  in  whitening  a  spectra, 
it  may  be  useful  for  detecting  radar  signals  that  arc  often 
obscured  by  powerful  but  anisotropic  radar  clutter.  The  MEM 
is  demonstrated  to  indeed  bo  a  likely  candidate  lor  further 
study  as  a  detector  of  weak  radar  signals  present  in  a  doppler 
radar  receiver. 

2 . 0  Radar  Data 

The  MEM  and  MJ,M  technique's  are  applied  to  actual 
radar  data  that  is  collected  by  a  14  sensor  linear  array 
within  a  laboratory  environment.  Two  si  ureas  ,  each  having 
power  levels  of  4  3  dB  above  the  averaged  background,  radiation, 
are  located  at  ly  and  22  degrees  from  the  0  degree  radial 
which  is  normal  to  th.  linear  array.  The  sensor  array  and 
radar  sources  are  both  located  within  an  anechoic  room  which 
substantially  reduces  specular  wall  reflections.  However, 
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the  non-specular  reflections  which  are  detected  by  the  array 
make  up  a  highly  colored  background  wavenumber  spectrum. 
Receiver  noise  and  sensor  phase  distortion  are  present,  bur 
are  unknown  spectral  quantities. 

The  14  sensor  array  data  is  collected  as  (I-Q  pair) 
data  so  that  14  complex  numbers  or  28  channels  of  data  are 
collectea  and  recorded.  The  radar  signals  are  short  duration 
pulses  so  that  each  set  of  14  complex  numbers  represent  one 
spatial  snapshot.  A  total  of  1024  snapshot  data  are  recorded. 

A  phased  array  antenna  pattern  has  a  mainlobe  beam- 
width  of  13.4  degrees  (at.  the  -3  dB  levels)  for  a  14  senscr, 
7.5  wavelength  array.  Consequently,  two  sources  separated  by 
only  4  degrees  are  not  resolved  with  a  conventional  summed, 
phased  array,  antenna  pattern. 

2 . 1  The  Maximum  Entropy  Method  (MEM ) 

Two  methods  of  obtaining  stable  MEM  wavenumber 
spectra  are  employed  for  analyzing  the  array  radar  data. 

One  method  employs  averaged  covariance  matricies  ^nd  the  other 
method  utilizes  averaged  prediction  error  filter  weights. 

Both  methods  provide  stable  MEM  wavenumber  spectra  with  suffi¬ 
cient  resolution  to  identify  the  two  closely  adjacent  sources. 
Both  averaging  methods  are.  described  and  demonstrated  in  an 
earlier  progress  report  dated  1,  Nov.  19  7')  (1)  .  The  MEM 

algorithm  employed  in  this  analysis  utilizes  the  Burg  tech¬ 
nique,  which  substantially  reduces  the  computing  effort 
required  oh  the  conventional  matrix  equation  formulation. 

The  Burg  technique  is  also  described  in  detail  in  the  earlier 
progress  report  (1) . 

2 .  2  The  Maximum  T.ikelihood  Method  (MLM) 

For  comparison  purposes,  the  MLM  is  also  used  to 
analyze  the  array  data.  The  MLM  is  evaluated  using  a  predic¬ 
tion  error  formulation  proposed  by  Burg  (2) ,  which  is  similar 
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to  the  MEM  spectral  formulation.  Because  of  the  similarity 
in  formulation,  the  same  prediction  error  filter  weights  ars 
computed  and  used  for  both  the  MEM  and  MLM  formulation.  The 
Burg  technique  is  also  employed  in  evaluating  the  MEM.  The 
MLM  formulation  is  as  follows: 


PN(kl  N  „  -l 


p  ‘  rr  (k) 

L  n  n 


n=l 


(1) 


where 


n+1  n 


1  -  (Yn+1)  (Yn+V 

1  '  n+r  '  'n  +  1' 


!  2) 


P.  =  r2  (autocorz'elation  of  data  samples) 

X  0 


r  (k)  =  Z  y n  eiKS(Ax) 
n  s-1  3 


(3) 


Ys  are  prediction  error  filter  weights 


k  =  (2tt/A)  sin  (0)  (wavenumber  component) 
8  =  signal  angle  of  incidence 

A  -■=  signal  wavelength 


2  .  3  Results  o  c_  Ana  lysis 


data 

digs 


The 
snapshot 
(1)  and 


MEM  and  MLM  wavenumber  spectra  for  the  first 
(first  set  of  14  complex  numbers)  are  shown  in 
(2)  respectively.  The  MEM  spectra  of  Fig.  (1) 
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is  typical  of  such  single  snapshot  spectra  in  that  it  contains 
many  very  sharp  peans  of  varying  magnitude.  In  such  spectra' 
it  is  often  difficult  to  determine  which  peaks  represent  signals. 
Hcwever,  there  are  three  large  peaks  in  the  MEM  snapshot  spectra 
that  are  at  azimuth  angles  in  the  vicinity  of  the  two  signal 
incidence  angles.  Ths  two  signals  are  located  at  angles  of 
18  and  2  2  decrees.  More  accurate  and  reliable  MEM  wavenumber 
spectra  are  obtained  if  particular  MEM  parameters  are  averaged 
over  several  sets  of  snapshot  data.  Subsequent  MEM  wavenumber 
spectra,  which  are  shown,  do  incorporate  averaging. 

The  MLM  snapshot  wavenumber  spectra  is  shown  in 
Fig.  (2)  where  one  large  and  prominent  peak  appears  at  an 
angle  in  the  vicinity  of  the  two  actual  signals.  In  only  one 
snapshot  of  data  the  MLM  spectra  is  observed  to  contain  very 
weak  extraneous  side  (noise)  peaks,  but  does  not  indicate  the 
presence  of  both  signals.  Averaging  is  incorporated  in  subse¬ 
quent  computed  MLM  spectra  with  the  hope  of  improving  resolu¬ 
tion  in  the  MLM  spectra. 

Other  MEM  and  MLM  wavenumber  spectra  are  shown  in 
Figs.  (3)  and  (4)  for  which  the  covariance  matrix  is  averaged 
over  6  sets  of  snapshot  data.  The  MEM  spectra  of  Fig.  (3)  does 
indicate  the  presence  of  two  signals  although  one  signai  peak 
is  about  10  or  12  dB  smaller  than  the  other.  The  side  peaks 
are  reduced  in  number  and  in  magnitude  in  comparison  with  the 
MEM  snapshot  of  Fig.  (1).  The  averaging  of  the  6  covariance 
matrices  has  improved  the  accuracy  and  reliability  of  the  MEM 
spectra,  but  the  resolution  remains  insufficient.  The  MLM 
spectra,  which  is  computed  for  the  covariance  matrix  averaged 
over  6  sets  of  snapshot  data,  has  not  changed  in  any  signifi¬ 
cant  way  in  comparison  with  the  snapshot  spectra  shown  in 
Fi  g.  (?)  , 

The  MEM  wavenumber  spectra,  which  is  shown  in  Pig, 

5,  is  computed  using  prediction  error  filter  weights  averaged 
ever  6  sets  of  snapshot  data.  The  MEM  spectra  of  Fig.  5  does 
indicate  the  presence  of  two  closely  adjacent  signals,  but  the 


¥ 

•V 


m 


-  '.u 1 B-  wnuupigj  uiiij 


resolution  capability  appears  inadequate  to  accurately  iden¬ 
tify  the  two  signals.  Side  peaks  are  still  quite  numerous 
and  large  in  magnitude  when  compared  with  the  snapshot  MEM 
spectra  of  Fig.  (1) . 

Averaging  does  improve  the  resolution  capability  of 
the  MEM.  The  MEM  spectra  is  shown  in  Fig.  (6)  where  the  pre¬ 
diction  error  filter  weights  are  averaged  over  12  sets  of 
snapshot  data.  In  Fig.  (6)  the  two  signals  are  well  resolved 
and  accurately  located  at  azimuths  of  18  and  22  degrees. 

Also,  side  peak  levels  are  reduced  but  are  not  suppressed 
as  well  as  with  the  MLM.  One  signal  peak  (at  18  degrees)  is 
about  1  dB  below  the  other  signal  peak  although  both  signals 
have  the  same  power  level. 

In  order  to  resolve  the  two  signals  with  MEM  spectra 
computed  with  averaged  covariance  matrices,  further  averaging 
is  necessary.  In  Fig.  (7)  the  two  signals  are  resolved  and 
accurately  located  by  averaging  the  covariance  matrix  of  an 
MEM  spectra  over  18  consecutive  sets  of  snapshot  data.  However, 
as  in  Fig.  (6) ,  one  signal  peak  is  about  1.5  dB  below  the 
other  signal  peak.  The  MLM  spectra  computed  using  the  same 
averaged  covariance  matrix  is  shown  in  Fig.  (8).  The  two 
signals  are  not  resolved  in  the  MLM  spectra,  and  are  not 
resolved  even  with  further  averaging  of  the  covariance  matrix 
in  which  27  sets  of  snapshot  data  are  utilized. 

3 . 0  Signal  Detection 

Because  of  the  spectral  whitening  capability  of 
the  maximum  entropy  method,  tnere  is  a  large  improvement  in 
the  ratio  of  SMR  (out)  to  SNR  (in) .  In  order  to  demonstrate 
the  detection  capability  of  the  MEM,  a  tvpical  signal- to¬ 
ol  uftei  environment  is  simulated  for  a  doppier  radar  receiver. 
Doppler  radar  clutter  is  simulated  as  shown  in  Fig.  (9)  where 
random  phased  clutter  bands  ha  ’  power  levels  typical  of  ground, 
rain,  and  interference  cluttei  Such  a  clutter  model  has  been 
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used  previously  by  Sawyers  (3)  in  his  demonstration  of  adapt¬ 
ive  filtering.  A  signal  having  a  0  dB  SNR  is  located  between 
the  clutter  at  the  frequency  ratio  of  . 375  as  denoted  by  the 
arrow  in  fig.  (9) .  The  signal  is  detected  as  shown  in  fig. 

(10)  by  applying  MEM  to  several  sets  of  32  data  samples  and 
using  24  filter  weights.  The  strong  clutter  bands  are  very 
effectively  whitened  by  MEM  such  that  the  largest  background 
peak  in  fig.  (10)  is  about  10  dB  below  the  signal  peak  level. 
Similar  results  ruay  be  obtained  for  any  signal  location. 

For  example,  in  fig.  (11)  a  signal  located  at  the  center  of 
the  interference  clutter  (.65)  is  equally  well  detected  again 
with  MESA  applied  to  consecutive  sets  of  32  data  samples 
using  26  filter  weights.  In  both  fig.  (10)  and  (11)  the  MEM 
filter  weights  are  averaged  over  30  consecutive  sets  of  32 
data  samples.  While  considerable  averaging  is  used  to  achieve 
the  results  indicated  in  figs.  (10)  and  (11)  ,  less  averaging 
of  fewer  filter  weights  may  also  achieve  satisfactory  signal 
detection,  but  with  less  resolution  capability. 

4 . 0  Summary 

The  maximum  entropy  method  is  effective  in  resolving 
signals  located  within  the  mainlobe  of  a  conventional  antenna 
pattern.  However,  the  MEM  is  accurate  and  stable  only  if  some 
form  of  parameter  averaging'  is  utilized.  As  demonstrated,  the 
MEM  as  compared  to  the  MLM  has  superior  resolution,  but  the 
MT.M  appears  to  more  effectively  whiten  the  background  (noise) 
spectra.  Previous  progress  reports  issued  under  this  investi¬ 
gation  contained  MEM  spectra  which  were  computed  using  only 
simulated  data.  In  this  report  all  MEM  spectra  are  obtained 
using  actual  radar  signals  received  by  a  14  sensor  linear  array 
Previous  simulated  data  contained  simulated  Gaussian,  white 
noise,  but  in  this  report  the  collected  radar  data  contains 
typical  receiver  and  sensor  noise  and  dj. sortion.  It  is  re¬ 
assuring  to  know  that  the  MEM  is  effective  in  resolving  closely 
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adjacent  signals  mixed  with  either  simulated  white  noise  or 
with  actual  and  typical  electronic  system  noise  spectra. 

The  results  achieved  with  the  MESA  doppler  radar 
detector  are  most  promising.  However,  such  exceptional 
results  did  require  an  average  of  over  32  sets  of  computed 
prediction  error  filter  weights.  Real  time  applications  may 
not  allow  for  such  a  large  average.  Decause  of  che  excellent 
spectral  whitening  characteristics  that  are  exhibited  by  the 
MLM,  perhaps  the  MLM  is  a  more  promising  technique  for  real 
tine  doppler  radar  detection.  Either  MESA  or  MLM  appear  to 
offer  considerable  improvement  in  radar  signal  detection. 
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Introduction 


Random  data  sampling  may  sometimes  be  a  useful  sampling 
technique  when  estimating  power  spectra  (1).  For  example,  the 
cost  of  a  large  antenna  array  may  be  significantly  .Lowered,  by 
reducing  the  number  of  antennas  (data  samples)  and  intentionally 
undersampling  through  the  use  of  a  random  distribution  of  anten¬ 
nas.  Random  (Lata  sampling  may  also  be  required  if  data  samples 
necessarily  have  a  random  distribution.  For  example,  radar  an¬ 
tennas  or  acoustical  sonobouys  may  be  dropped  from  the  air  and 
permitted  to  fall  freely  to  the  ground  or  water,  and  land  in  a 
random  planar  distribution. 

With  use  of  random  data  sampling,  the  side  lobes  of  a  Four¬ 
ier  power  spectra  have  a  random  distribution  in  both  number  and 
location.  However,  thi  typical  Fourier  spectral  characteristics 
such  as  resolution  and  side  lobe  levels  may  be  retained  in  an  av¬ 
erage  of  many  such  Fourier  power  spectra  computed  using  random 
sampling.  Useful  Fourier  power  spectra  also  result  from  random 
data  sampling  even  with  very  sparse  (undersampled)  data  sots. 

It  is  questionable  whether  random  data  sampling  is  effec¬ 
tive  when  the  maximum  entropy  method  (MEM)  is  used  to  compute  an 
estimated  power  spectra.  Examples  of  both  the  Four  i.ei  and  MEM 
power  spectra  are  computed  using  random  data  sampling  to  demon¬ 
strate  the  effectiveness  of  random  data  sampling  with  these  two 
power  spectra  estimation  techniques , 

Examples  of  Fourier  A ntenna  Patterns 


Fourier  antenna  patterns  may  be  computed  by  taking  the 
spatial  Fourier  transform  c£  signals  incident  to  an  "  itenna  array. 
The  Fourier  transform  of  N  spatial  data  points  y  is  as  follows: 


X  (k) 
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n-  1 


-i  kx 
e  n 
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where  y  is  the  measured  data  and  k  Ls  the  signal  wavenumber  de¬ 
fined  as  follows: 


k  =  (2tt/A)  s.i.n(Q) 


*  ik„x 
y  =  A  e  On 
'  n 


t  n  e 


iij) 


for 


k  Q  =  C’ir/A)  sin  ( 0  0  ) 


'I’he  amplitude,  wavelength,  and  angle  of  the  incident  signal  in 
respectively  given  by  A,  A,  and  0.  The  noise  amplitude  g  has  a 
Gaussian  distribution,  and  the  noise  phase  '(>  has  a  wh.i  to  distri¬ 
bution.  The  N  array  antennas  are  distributed  at  positions  x  along 
a  one  d.i.me.ns  i.ona.I.  axis.  Tn  a  random  l.y  d:i  sl.r  ibuted  antenna  array 
the  positions  x  are  located  with  equal  probab.1.1  i  ty  along  a  lone 
of  length  I.  according  to  a  white  distri.Jjutj.en  function. 

The  resulting  Fourier  antenna  patterns,  computed  using  ran¬ 


dom  sampling  are  shown  in  Figs.  1.-3.  In  the  fit  rd  example  ol  Fig. 

1  a  signal  is  detected  with  1 6  randomly  distributed  antennas  loca¬ 
ted  in  a  straight  line  .array  having  a  total  length  of  H  wave  lengths . 
Side  lobe  structure  is  not.icubl.y  higher  than  those  obtained  with 
hall: -wavelength  sumpl  ng.  However,  in  an  average  of  Ml  such  an¬ 
tenna  patterns,  the  result  in'  averaged  pattern,  shown,  in  F.iu.  2, 
lias  reduced  side  lobe  levels  that  are  about  l(.)  dll  below  the  main 
l.ohe.  Further  averaging  of  such  antenna,  patterns  will  result  in 
lower  side  l.obe  levels  approaching  the  -13.6  dB  minimum  level. 

The  width  of  the  main  lobe  in  examples  shown  in  Figs.  .1.  and 


sampling.  Only  the  antenna  pattern  side  lobe  structure  is  sign.i.1  i- 
euntly  altered  when  using  random  data  sampling .  In  a  third  exam¬ 
ple  an  average  of  10  antenna  patterns  is  shown  in  Fig.  3,  where 
each  Fourier  antenna  pattern  is  computed  using  a  very  sparse  data 
set  cons  Luting  of  only  4  data  samples  ,  which  are  randomly  distrib¬ 
uted  over  a  single  dimension  of  0  wavelengths  in  length. 


The  re- 
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sultant  averaged  antenna  pattern  is  surprisingly  good,  v.hen  it  is 
considered  that  the  detected  signal  is  unde. r sampled  by  a  factor  of 
four . 

Examples  of  MEM  Wavenumber  Patterns 

Because  the  maximum  entropy  method  may  he  derived  using 
the  discrete  Wiener  prediction  filter  (2),  it  is  unlikely  that 
meaningful  data  predictions  are  possible  with  random  data  sampl¬ 
ing.  For  example  a  predicted  signal  ft  is  given  by  the  discrete 
convolution  as  follows: 


N 

ft  =  >:  b  x 

m  .  n  m-n 

nr.  1 


CD 


where  the  data  set  x  must  be  a  set  of  uniformly  sampled  data. 

m 

If  only  a  relatively  small  number  of  data  samples  are  massing,  then 
oqn.  3  still  may  provide  useful  predicted  data  .from  which  the  pre¬ 
diction  error  may  be  minimized.  However,  it  the  entire:  data  set 
is  obtained  using  random  sampling,  then  the  data  net  ft  predicted 
using  oqn.  3,  is  not  expected  to  be  very  meaningful.  Hineo  the 

prediction  f  i.  Iter  b  is  derived  from  a  minimum  error  derived  with 

n 

eqn.  3,  the  prediction  filter  b  may  also  be  an  inefficient',  filter 

if  the  data  set  x  is  derived  using  random  data  sampling. 

The  MEM  all -pole  wavenumber  spectra  formula! ion  does  not 

present,  any  apparent  problem  with  use  of  random  data  sampling,  since 

the  MEM  wavenumber  spectra  I.  or  mu'!  a  t  ion  eon  ha  i  no  only  a  discrete 

Fourier  transform  of  the  prediction  filter  b  as  indicated: 

L  n 

M  _  .  ,  2 

F(k)  -  (P  /k)  /  |  I  “  >3  b  e  '  "  Xfi  |  (4) 

n  1 


However,  since  the  prediction  f  il  ter  b  is  .likely  to  be  in  error  , 
the  MEM  power  spectra  as  given  by  eqn.  1  may  also  not  be  very 
meanin  [fill  when  computed  using  random  data  sampling.  Example.;  u( 
MEM  wavenumber  spect  ra  are  demonstrated  in  Fi  gs.  4  -f» ,  where  the 
first  example  is  a  MEM  snapshot  wavenumber  spectra  computed  with 
uniform  data  sampling.  Other  MEM  wavenumber  spectra  shown  arc 


1 


are  computed  using  random  data  sampling. 

The  signal  is  accurately  detected  in  Fig.  4  by  an  8  wave¬ 
length  linear  antenna  array  having  uniform  half-wavelength  antenna 
spacing.  The  MEM  wavenumber  spectra  shown  in  Fig.  5  is  computed 
using  16  antennas  randomly  spaced  along  a  straight  line  that  is  8 
wavelengths  long.  While  both  snapshot  wavenumber  spectra  shown  in 
Figs.  4  and  5  are  computed  for  a  signal  incident  to  the  antenna 
array  at  30  degrees  and  with  SNR  of  20  dB,  only  the  MEM  spectra 
computed  with  uniform  half-wavelength  sampling  accurately  detects 
the  signal.  The  MEM  wavenumber  spectra  shown  in  Fig.  5  indicates 
that  a  signal  is  detected  at  about  16  degrees.  Other  MEM  spectra 
computed  with  random  data  sampling  (not  shown)  als<  contained 
similar  false  alarms.  An  average  of  6  such  MEM  wavenumber  spectra, 
where  each  spectra  is  computed  using  random  data  sampling,  is 
shown  i  i  Fig.  6.  Even  though  16  antennas  are  located  randomly 
along  a  straight  line  having  a  length  of  8  wavelengths,  neither 
the  individual  6  computed  spectra  or  the  average  spectra  indicated 
a  signal  detection  in  a  6  degree  window  about  the  signal  angle  of 
30  degrees.  The  averaged  MEM  wavenumber  spectra  shown  in  Fig.  6 
contains  only  false  alarm  peaks. 

As  anticipated,  the  MEM  wavenumber  spectra  computed  using 
random  data  sampling  are  not  useful  spectra,  since  MEM  does  require 
the  use  of  uniform  data  sampling  at  minimal  half -wave length  inter¬ 
vals.  Of  course  if  it  is  possible  to  obtain  accurate  estimates 
of  a  uniformly  spaced  data  set  using  some  extrapolation  method- 
then  useful  MEM  spectra  could  conceivably  be  computed  from  a  ran¬ 
dom  data  set.  It  is  doubtful  that  any  useful  spectral  estimates 
can  be  obtained  with  any  spectral  estimation  technique  based  upon 
the  prediction  filter,  if  random  data  sampling  is  utilized.  How¬ 
ever  conventional  Fourier  spectral  estimation  techniques  do  provide 
useful  spectra  using  randomly  sampled  data. 
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RANDOMLY  LOCATED  ANTENNAS 


THE  W-K  SPECTRAL  ESTIMATION  TECHNIQUE 


Introduction 

The  maximum  entropy  spectral  analysis  (MESA)  tech¬ 
nique  and  the  autoregressive  (AR)  spectral  analysis  technique 
are  both  limited  in  accuracy  and  stability  cy  increasingly 
noisy  data.  The  MESA  and  AR  techniques  may  be  unnecessarily 
noise  limited  due  to  the  definition  of  the  Wiener  prediction 
error,  which  is  minimized  in  both  techniques.  It  has  been  noted 
(1)  that  the  conventional  definition  of  the  prediction  error 
actually  contains  two  error  components.  One  error  component 
is  the  actual  prediction  error  inherent  in  the  predicted  data 
set  xm  consisting  of  M  data  points.  The  other  component  error 

is  the  noise  present  in  the  actual  M  data  samples  x  . 

r  m 

The  prediction  error  may  be  redefined  so  that  there 
is  only  one  error  component,  i.e.  the  inherent  prediction  error. 

By  redefining  the  prediction  error,  it  is  anticipated  that,  spec¬ 
tral  estimation  accuracy  and  stability  of  the  resulting  spectral 
estimator  will  be  substantially  improved  at  the  lower  signal-to- 
noi se  power  levels.  The  resultant,  new,  spectral  estimation  tech¬ 
nique  is  referred  to  as  the  Wiener-King  (W-K)  spectral  estimator. 

Prediction  Error  Definition 

The  conventional  definition  for  the  time  dependent 

prediction  error  e  is  as  follows: 

1  m 

e  =  x  -  x  for  (m=l,M)  (1) 

in  m  in 

where  there  are  M  data  samples  x  and  M  predicted  data  points 

m 

x  .  The  predicted  data  points  are  defined  as  follows: 
m  ^ 

N 

x  =  T,  a  x  for  N  <  M  (2) 

m  .  n  m-n 

n=l 

where  there  are  N  prediction  filter  coefficients  a^  (n=l,N). 
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Thp  conventional  prediction  error  e 

rr 

of  the  noisv  data  as  well  as  the  predicted  data.  The  data  noise 
may  be  readily  eliminated  by  redefining  the  prediction  error  as 
a  function  of  noise  free  spectral  components.  A  spectral  com¬ 
ponent  yP ,  which  has  frequency  f  is  defined  by  the  exponential 
function  as  follows: 


p  i2fTf  m(At) 
yl  =  e  P 
■ 1  m 


The  "new"  noise  free  prediction  error  e1^  is  defined  as  follows: 


x 


m 


eP  =  ei2irfpmUt) 


N 

Z  a 
n=l 


nxm-n 


(4) 


The  "new"  prediction  error  contains  only  one  error  component, 
which  is  the  inherent  error  present  in  the  prediction  function 
x  .  A  total  mean  squared  prediction  error  may  be  computed  by 
summing  the  time  depend-. .L  squared  error  computed  for  all  M  pre¬ 
dicted  data  points. 


Prediction  Filter 


A  "new" 
mi  zing  a  redefined 
total  mean  squared 


prediction  filter  may  be  defined  by  mini¬ 
total  mean  squared  prediction  error.  A 
prediction  error  Ep  may  be.  *  fined  as  follows  : 


o 

e" 

(ep)*  ^ 

=  (VP 

x  ) 

* 

m 

m 

•*  m 

m 

*  m 

m 

■  vS< 

vP) 

*  ^ 

:<  ( 
rn 

y  )  *  - 
J  m 

X*vP 
\\\x  m 

-P 

-  (1/ 

(  2M+  1 

)  ) 

V 

( c 

P), 

.  m 

-  n  7  — 


The  total  mean  squared  prediction  error  is  expressed  in  detail 
by  combining  eqns.  (5}  and. (6)  as  follows; 


EP  = 


1/ ( 2M+ 1 ) 


M  d  M 

E  vP(yP)*  -  E  a  E  x  (yp)* 
-m  ■'m  ,  n  ..  m-n  m 

m=-M  n=l  m=-M 


N  M  N  m  M  i 

E  a*  E  x*  yp  +  E  Z  a  a*  E  x  x*  I 

,  n  ..  m-n  m  ,  ,  ,  n  k  m-n  m-k 

n=l  m=-M  k=l  n=l  m=-M 


(7) 


By  utilizing  the  definition  of  the  autocorrelation  and  cross- 
correlation  functions,  expression  (7)  is  simplified  as  follows: 


E 


P 

N 


M 

E  a 
n=l 


rp 

n  xy 


(n) 


N  N 
+  T  E  a 
k=l  n«l 


a*r 
n  k  x 


(k-n) 


(8) 


where  the  autocorrelation 
tion  of  function  x  is  r  . , 
tions  yP  and  x  is  r  and 
The  total  mean 
mized  with  resoect  to  the 


of  function  yp  is  rp,  the  autocorrela- 
and  the  crosscorrelation  of  the  func- 
r 

yx 

squared  prediction  error  may  be  mini¬ 
unknown  filter  coefficients  a  as  follows: 

n 


r h  (:.)  + 
yx 


N 

E  a  r  (n-i )  =  0 

n  a 


(9) 


v/here  1  =  i ,  N . 


The 

tions  that  may 
cients.  The  N 


resultant  expression  eqn.  (9)  is  a  set  of  N  equa¬ 
te  simultaneously  solved  for  the  M  filter  coeffi- 
equations  may  be  written  as  follows: 
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1  =  1, 

alrxc0) 

+  a2rx(1) 

+  ‘ *  *  +  aNrx (v-1) 

=  rp  (1) 
yx 

r\i 

II 

o? 

Vx1'11 

+  a  2  r  ,t  ( 0 ) 

+  1  aMrx  (>’~2 ) 

=  rpy(2) 

=  N,  air.<  (-N+1)  +  a2r  (-N+2)  + 


+  aNrx(0) 


=  rp  (N) 


The  N  equations  that  define  the  prediction  filter  coefficients 
an  may  also  be  expressed  in  matrix  notation  as  follows: 


rxC0)  rx(l) 

' •  r  (K-l) 

t - 

1— 1 

X 

0-i  >1 

ih _ 

r x  ( 1 )  rx(0i 

■  • 

•  . 

*’*  rx(V-2) 

rp  (2) 
yx 

r  ( N- 1 )  r.  (N-2) 

L 

I 

O 

X 

u 

rP  (N) 
yx 

where  it  is  noted  that  r  ( —  n )  =  r  (n) 

X  X 


The  matrix  equation  gi”en  by  eqn.  (11)  may  be 
written  in  shorthand  notation  as  follows: 


a  r  =  r*- 
x  vx 


The  solution  for  the  set  of  filter  coefficients  is  obtained 
by  inverting  the  autocorrelation  matrix  r  and  solving  for 
the  coefficients  as  follows : 


(13) 


It  is  necessary  to  solve  eqn.  (13)  for  all  coefficients 

(a.  ,  a0,  a,, )  for  everv  spectral  comDonent  f  of  inter- 

12  N  £  p 

est.  However  it  is  possible  to  determine  spectral  regions 
of  interest  by  simply  computing  the  crosscorrelation  vec¬ 
tors  at  every  frequency  f  .  In  this  way  it  is  necessary 

only  to  solve  for  the  prediction  filter  coefficients  in 
spectral  regions  where  the  crosscorrelation  vectors  indicate 
a  strong  correlation. 

Power  Spectra 

Consider  another  prediction  error  e  which  is  de- 

m 

fined  to  include  the  unknown  signal  amplitude  A  ,  i.e. 

P 


r  = 

A 

ep 

m 

P 

m 

„  __ 

A 

yp  - 

X  I 

P 

mj 

r  m 

M 

£  ~ 

A 

|yp  - 

Z  a 

m 

P 

l> 

n=l 

transform 

of 

eqn . 

(15) 

for  the  spectral  amplitude  A_^  as  follows: 


f(f)  =  A  5  ( f  -  f  )  -  X(f)  l  a  el2''Tfn(At' 

p  P  n=l  n 


,  v  t  £  \  r  i2u  n  ht 
1  -  X  ( f„  )  E  a  e  p 
P  .  n 
n=l 


The  power  spectra  P(f  )  is  given  by  the  following 


expression : 


-  loo- 


2 


(17) 


P(f  ) 
P 


L  -  X(fp)  j^n6 


2irf (At) 


where  5(f),  which  is  the  Fourier  transform  of  the  prediction 
error  e  ,  may  be  evaluated  in  the  spectral  region  void  of  sig¬ 
nals,  i.e.  £  /  £  .  For  the  present  it  useful  to  consider  the 
P 

all  pole  representation  of  the  W-K  power  spectral  technique  by 

evaluating  the  ratio  of  P(f  ) /%^ (£  ). 

P  P 

It  is  anticipated  that  the  W-K  all  pole  power  ratio  spectra 

computed  with  eqn.  (17)  will  be  more  stable  but  just  as  accurate 

as  comparative  spectra  computed  with  MESA.  The  formulation  of 

the  W-K  spectral  estimator  as  given  by  eqn.  (17)  differs  from 

that  of  MESA  in  that  the  Fourier  transform  X(fp)  of  the  data  set 

is  utilized  as  indicated  in  eqn.  (17) .  Of  course  the  prediction 

filter  coefficients  a  are  different  from  the  prediction  coef- 

n 

ficients  used  in  MESA. 

The  presence  of  the  data  transfer  function  X(f  )  is  a  sta- 

p 

blizing  factor  in  the  power  spectral  eqn.  (17).  Spectral  reso¬ 
lution  and  accuracy  of  the  W-K  spectral  estimator  are  expected 
to  be  comparable  to  similar  characteristics  of  MESA  r  since  both 
snectral  estimation  methods  have  a  similar  all  pole  formuliza- 
tion.  However  the  W-K  spectra],  estimator  is  expected  to  have 
better  stability  and  whiter  spectra  than  that  obtainable  with 
use  of  the  MESA  tecnnique. 

Example  Wavenumber  Power  Spectra 

A  W— K  wa veil  Umber  power  spectra  may  be  obtained  Crum  eqn. 
(17)  by  introducing  the  wavenumber  k  and  spatial  dimension  x 
with  aid  of  the  variable  relationships  that  follow: 

k  =  V  (  2  t  f  ; 
x  -  t 

where  V  is  the  wave  velocity.  With  the  indicated  change  of 
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variables,  examples  of  W-K  wavenumber  power  spectra  are  compu¬ 
ted  for  only  a  few  simple  signal  and  noise  conditions.  At  this 
time  the  W-K  spectral  estimation  technique  is  evaluated  only  to 
indicate  that  the  W-K  spectral  estimator  is  a  viable  and  useful 
method  for  obtaining  stable,  high  resolution  spectral  estimations 
with  only  a  single  (snapshot)  set  of  data. 

A  single  signal  is  shown  detected  at  28.5  degrees  in  Fig.  1 
using  a  181  point  W-K  wavenumber  power  (ratio)  spectra  as  given 
by  eqn.  (17) .  The  simulated  signal  is  actually  incident  to  an  8 
element  linear  antenna  array  at  an  angle  of  30  degrees  and  with 
an  SNR  of  20  dB.  The  W-K  spectrum  shown  in  Fig.  1  is  computed 
usinn  only  3  filter  coefficients.  The  snapshot  spectra  is  well 
whitened  by  the  W-K  estimation  technique  although  the  signal 
location  is  imprecise  with  an  error  of  1.5  degrees.  The  same  sig¬ 
nal  and  noise  model  is  used  to  obtain  another  W-K  wavenumber 
spectra  using  4  filter  coefficients  as  shown  in  Fig.  2.  Again 
the  W-K  spectra  is  well  whitened,  and  again  the  signal  peak  is 
located  with  the  same  imprecision.  A  comparison  of  these  two 
initial  W-K  spectra  indicate  a  smaller  filter  size  may  provide 
better  noise  suppression  by  almost  10  dB  in  the  examples  given. 

For  comparison  a  MESA  wavenumber  spectra  is  computed  for 
the.  same  sign.'  1  and  noise  model  used  in  the  two  W-K  spectra  ex¬ 
amples.  However  five  filter  coefficients  are  used  in  the  compu¬ 
ted  MESA  example  to  better  i] lustrate  the  large  noise  peaks 
that  frequently  appear  in  MESA  snapshot  spectra.  It  should  be 
noted  that  whitened  spectra  do  result  when  MESA  is  applied  to 
covariance  matrix  data  that  has  been  averaged  over  several  snap¬ 
shots  of  data.  The  example  MESA  wavenumber  spectra  is  shown  in 
Fig.  3,  where  indeed  5  large  peaks  appear  in  the  spectrum,  but 
only  the  peak  located  at  about  31,5  degrees  is  representative 
of  the  signal. 

In  another  example  two  signals  incident  to  the  antenna 
array  at  -30  and  30  degrees  are  shown  correctly  detected  by  the 
W-K  spectral  estimator  in  Fig.  4.  The  overall  W-K  spectra  is 
well  whitened  with  use  of  5  filter  coefficients.  For  comparison 
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SIGNALS  AT  -30  and  30  DEGREES 


the  same  two  signals  (SNR  =  20  dB)  and  noise  are  shown  in  Fig. 

5  detected  by  a  conventional  Fourier  wavenumber  antenna  pattern. 
The  strong  Fourier  side  lobe  structure  that  is  present  in  Fig.  5 
is  only  very  mildly  apparent  in  the  W-K  all-pole  wavenumber  spec¬ 
tra  shown  in  Fig.  4.  If  the  two  spectra  of  Figs.  4  ana  5  are 
superimposed,  the  Fourier  side  lobe  structure  is  found  to  occur 
at  the  same  spectral  locations  of  the  W-K  side  peaks.  This  simi¬ 
larity  in  side  peak  location  is  an  indication  of  the  stablizing 
influence  brought  about  by  the  presence  of  the  Fourier  transfer 
function  X(f  )  that  appears  in  the  W-K  spectra  formalization. 

The  resolution  capability  of  the  W-K  spectral  estimator 
is  illustrated  in  Fig.  6,  where  two  closely  adjacent  signals 
(6  degrees  apart)  are  resolved  by  the  W-K  spectral  estimator. 

For  comparison  purposes,  the  same  two  signals  (SNR  ~  20  dB)  are 
shown  detected  in  Fig.  t  as  a  single  peak  in  a  conventional 
Fourier  antenna  pattern. 

As  had  been  anticipated,  the  W-K  all  pole  spectral  esti¬ 
mator  is  a  very  stable ,,  high  resolution  spectral  estimation 
technique.  In  all  examples  shown  the  computed  wavenumber  spectra 
is  well  whitened  by  the  W-K  spectral  estimator.  Spectral 

whitening  is  always  achieved  with  only  one  snapshot  of  data. 
However  some  apparent  improvements  in  this  original  version  of 
the  W-K  spectral  estimation  technique  are  indicated  by  the  ex¬ 
ample  spectra.  The  small,  numerous  peaks  that  are  often  evi¬ 
dent  in  the  W-K  spectra  are  a  problem  when  they  occur  in  the  vi¬ 
cinity  of  a  signal  peak.  In  such  instances  the  presence  of  the 
small,  noisy  peaks  may  be  misintrepr eted  as  a  collection  of 
several  very  close  signals.  Perhaps  an  improved  version  of  the 
W-K  estimator,  or  a  simple  correction  of  some  overlooked  program¬ 
ming  error  may  serve  to  eliminate  or  at  least  minimise  the  pres¬ 
ence  of  the  small  annoy inq  noisy  peaks. 

The  derivation  of  the  W-K  spectral  estimator  has  util¬ 
ized  only  the  forward  prediction  error.  Perhaps  a  minimiza¬ 
tion  of  the  total  (forward  and  backward)  prediction  error  will 
serve  to  further  improve  the  W-K  SDectral  estimator. 
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Introduction 

The  Wiener  prediction  filter  is  useful  in  deriving 
the  high  resolution  maximum  entropy  and  W-K  spectral  estima¬ 
tion  techniques.  Both  of  these  spectral  estimators  employ 
prediction  filters  that  are  defined  by  minimizing  a  mean 
squared  Wiener  prediction  error.  The  evaluation  of  these 
prediction  filters  may  require  the  computation  and  manip¬ 
ulation  of  large  matricies ,  which  usually  requires  consid¬ 
erable  computational  time. 

While  these  spectral  estimation  techniques  have  very 
unique  and  useful  properties,  it  is  also  possible  to  obtain 
another  useful  spectral  estimator  by  utilizing  a  very  simple, 
and  quickly  computed,  prediction  filter.  The  resultant  spec¬ 
tral  estimator  is  named  the  "Fourier  spectral  estimator",  be¬ 
cause  it  employs  only  the  Fourier  transform  of  the  data  set. 
The  Fourier  spectral  estimator  (FSE)  may  be  derived  either  as 
an  all-pole  or  as  a  zero-pole  model.  Both  models  are  pre¬ 
sented  in  the  analysis  that  follows. 

Wiener  Prediction  Filter 

Narrowband  spectral  components  in  a  detected  signal 
may  he  characterized  with  use  of  the  Wiener  prediction  filter 
as  follows: 


A  c 
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iw  t 
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k(t) 


(t) 


(i) 


* 

#1 

■p 


where  A^  and  are  the  amplitude  and  frequency  of  a  narrow- 
band  component,  and  e(t)  is  the  spectral  prediction  error. 

The  prediction  signal  x(t)  may  be  taken  as  the  convolution 
of  the  prediction  filter  b(t)  with  the  measured  function  x(t) 
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as  follows: 


CO 

x(t)  =  /^b (T-t) X (t) dr 


(2) 


Filter  Definition 

The  prediction  filter  b(t)  is  usually  defined  by  mini¬ 
mizing  the  mean  squared  prediction  error.  However  it  is  much 
easier  to  choose  a  prediction  filter  that  yields  a  zero  pre¬ 
diction  error  at  the  exponential  frequency  u)  .  The  filter 
is  defined  in  the  frequency  domain  by  combining  the  Fourier 
transforms  of  eqns.  (1)  ari  (2)  as  follows: 

A  6  (oj~o)  )  =  B  (a>  ,t)  X  (co)  -  E  (co)  (3) 

IT  P 

If  E(co)-*0  as  ,  then  it  follows  from  eqn.  (3)  that 

P 

B  (oj  ,  t)  -  A  /  X  (co  )  (4) 

P  P  P 

All  Pole  Solution 

Even  though  eqn.  (4)  is  correct,  cop  is  not  known  and 
it  is  not  possible  to  evaluate  X(o  ).  However,  for  high  SNR 
values  X(o)p)  is  approximately  equal  to  its  maximum  magnitude  XQ, 
which  is  more  readily  determined.  The  prediction  filter  trans¬ 
fer  function  may  be  specified  approximately  as  follows  : 

B(Vt)  "  Ap  /  X0  (5) 

By  combining  eqns.  (3)  and  (5)  the  narrowband  component 
amplitude  is  given  as  follows: 

A?  *  Ap  X(cop)  /  XQ  -  E-(cop)  (6) 

Ap  *  -EMcop)  /  (1  -  X(«p)/X0}  P> 
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where  E' (co^)  is  the  prediction  error  due  to  the  approximate 
value  of  the  prediction  filter  transfer  function  used  in  eqn. 
(7),  An  all  pole,  white,  spectral  power  ratio  may  be  obtained 
from  eqn.  (7)  as  follows: 


|A  /E' |  =  1  /  jl  -  X(op)/XQ 


2 


(8) 


The  frequency  «  of  the  exponential  components  may  be  approxi¬ 
mated  by  determining  the  complex  poles  of  eqn.  (8) . 

Zero-Pole  Solution 

At  spectral  regions  other  than  «  eqn.  (3)  reduces  to 

P 

the  following  expression  for  the  prediction  error: 


E(oj)  =  B  (o)ft)  X(co) 


(9) 


The  most  obvious  choice  for  the  filter  transfer  function  is 
zero,  however  it  may  be  more  preferable  to  choose  B(co,t)  =  1.0 
so  that 


E  (to)  =  X(co)  (10) 

The  choice  of  eqn.  (10)  permits  other  filter  methods  to  be 
applied  for  the  elimination  of  strong  interference  signals, 
that  may  otherwise  appear  as  a  pole  in  eqn,  (8).  A  zero-pole 
Fourier  power  spectra  results  from  combining  eqn.  (10)  and 
eqn.  (7)  to  obtain  an  expression  useful  over  the  entire  fre¬ 
quency  domain.  The  zero-pole  Fourier  spectral  estimator  is 
qiven  as  follows: 

ApAp  ~  |X(to)  /  (1  X(o])/X0)|2  (11) 
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Example  Spectra 


The  conventional  Fourier  antenna  pattern  is  compared 
with  several  examples  of  the  stable,  all-pole  and  zero-pcle 
Fourier  spectral  estimator  (FSE) .  All  examples  are  snap¬ 
shot  antenna  and  wavenumber  spectral  patterns,  that  are  com¬ 
puted  for  signals  of  20  dB  SNR  which  are  incident  to  an  8 
element,  linear,  antenna  array.  The  white,  Gaussian  ampli¬ 
tude  noise  data  is  identical  (IR=3)  in  every  example  pre¬ 
sented. 

The  conventional  Fourier  antenna  pattern  for  one 
signal  incident  at  30  degrees  is  shown  in  Fig.  1.  The  same 
signal  is  shown  detected  with  the  all-pole,  Fourier  esti¬ 
mator  in  Fig.  2,  and  with  the  zero-pole  Fourier  estimator  in 
Fig.  3.  Other  comparative  examples  include  two  widely  spaced 
signals  incident  at  -30  and  30  degrees  that  are  shown  in  Figs. 
4,5,  and  6;  two  closely  adjacent  signals  incident  at  30  and 
45  degrees,  which  are  shown  in  Figs.  7,8,  and  9.  One  other 
set  of  examples  (Figs.  10,11,  and  12)  show  five  gnals,  all 
with  the  same  power  level  (20  dB  SNR) ,  incident  at  angles  of 
-60,  -45,  0  t  30,  and  37  degrees. 

In  all  the  examples  shown,  the  FSE  technique  appears 
to  have  improved  spectral  accuracy  and  improved  resolution  and 
detection  capability  in  comparison  with  the  conventional  Fourier 
power  spectral  method.  The  derivation  of  the  FSE  technique 
appears  to  indicate  that  these  improved  spectral  character¬ 
istics  are  more  significant  with  increasing  SNR  levels. 

In  the  Fourier  antenna  pattern  of  Fig.  10,  the  signal 
power  levels  of  the  three  signals  detected  near  0,  30,  and  37 
degrees  c  re  in  error  by  several  dB.  These  same  errors  are, 
in  effect,  magnified  by  the  FSE  examples  giver,  in  Figs.  11 
and  12.  Consequently  signal  power  levels  indicated  by  either 
F.'E  technique,  are  in  much  greater  error  than  the  power  levels 
indicated  by  the  conventional  Fourier  antenna  pattern.  It 
may  be  difficult  to  measure  either  absolute  or  relative  power 
levels  from  computed  FSF  wavenumber  spectra. 
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Perhaps  all  signal  parameters  may  be  more  accurately 
determined  by  using  both  the  corvrentional  Fourier  spectral 
and  the  Fourier  spectral  estimator  techniques.  Hopefully 
the  FSE  detection  ,  resolution,  and  spectral  characteristics 
may  be  determined  in  a  future  extensive  analysis. 
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